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ABSTRACT 


— fe.  LOWTRAN  and  MODTRAN  v/ero  evaluated  in  the  2. 0-5. 5  micron 
region  against  field  collection  data  at  high  zenith  angle/long 
path  lengths  to  determine  the  degree  of  uncertainty  associated 
with  these  models  under  these  conditions.  Matching  data  sets 
were  developed  using  data  from  the  Air  Force  Geophysics 
Laboratory  Flying  Infrared  Signatures  Technology  Aircraft 
(FISTA)  as  the  field  reference.  The  Kolmogorov-Smirnov  test 
was  applied  to  determine  the  degree  to  which  the  outputs  of 
LOWTRAN  and  MODTRAN  follow  the  same  distribution  as  the  field 
data.  The  percent  difference  between  the  model  and  field  data 
was  also  studied. 

Agreement  between  the  model  and  field  data  was  found  to 
be  better  than  97%  for  most  cases.  Median  percent  difference 
was  within  10%  for  zenith  angles  less  than  90  degrees . 
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1.  INTRODUCTION 


Proper  interpretation  of  events  viewed  through  the 
earth's  atmosphere  requires  that  the  effects  of  atmospheric 
transmission  and  radiance  be  quantified  when  analysis  of  the 
remotely  sensed  data  is  performed.  Computer  models  of  the 
atmosphere  have  been  developed  to  allow  system  designers  and 
data  analysts  to  perform  this  quantification.  Two  of  the  more 
commonly  used  models  are  LOWTRAN  and  MODTRAN,  developed  by  the 
USAF  Geophysics  Laboratory.  These  models  are  computer  codes 
which  calculate  atmospheric  transmission  and  background 
radiance  in  the  ultraviolet  (UV)  through  infrared  (IR)  regions 
of  the  electromagnetic  spectrum. 

Traditionally,  remote  sensing  has  been  performed  within 
a  viewing  cone  ±  45°  from  earth  normal,  either  uplooking  or 
downlooking.  This  self-imposed  limitation  by  remote  sensing 
system  designers  and  data  analysts  was  due  to  concerns  over 
optical  system  quality  and  atmospheric  effects.  The  increased 
call  for  remote  sensing  to  fill  needs  in  both  commercial 
(industrial  and  environmental)  and  military  areas  is  resulting 
in  systems  which  are  required  to  make  use  of  wider  viewing. 
Examples  of  this  are  the  desire  to  use  data  from 
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meteorological  satellites  for  use  in  nonmeteorological 
applications  (eg  like  LANDSAT  images)  and  broader  area 
coverage  from  environmental  missions  of  high  altitude  aircraft 
such  as  the  NASA  U-2. 

This  thesis  reviews  the  modeling  by  the  LOWTRAIJ  and 
MODTRAN  codes  and  compares  their  statistically  against  actual 
field  measurements  in  high  zenith  angle  viewing  (60-90 
degrees)  with  long  path  distances  (>50  km) .  By  testing  the 
codes  along  these  extreme  viewing  conditions  it  is  hoped  that 
the  effectiveness  of  the  modeling  can  be  ascertained.  Key  to 
this  focus  on  modeling  errors  is  the  use  of  high  altitude 
measurements.  By  using  high  altitude  measurements  atmospheric 
variability  is  significantly  reduced  (the  atmosphere  is 
relatively  stable  and  predictable  at  tropospheric  and 
stratospheric  altitudes) .  In  a  statistical  sense  this  allows 
the  effects  of  atmospheric  uncertainty  in  any  given 
measurement  to  be  blocked  out  and  only  modeling  errors  should 
remain. 

Determining  the  modeling  errors  is  important  for  two 
general  classes  of  problems.  In  many  remote  sensing 
applications  either  LOWTRAN  or  MODTRAN  are  used  to  determine 
the  effect  of  the  atmosphere  on  data  collected  in  a 
measurement  program.  This  effect  is  then  used  to  back  out  the 
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losses  to  arrive  at  the  original  radiance  of  the  object  of 
interest.  Modeling  errors  would  thus  lead  directly  to 
measurement  errors  in  the  processed  data.  These  errors  could 
lead  to  incorrect  results  in  the  analysis  of  that  data. 

The  codes  are  also  useful  in  this  problem's  inverse. 
Sensors  are  usually  designed  against  a  minimum  signal  of 
interest.  This  signal  is  often  defined  as  the  signal  of  the 
target  less  the  attenuation  of  the  atmosphere.  The 
attenuation  is  determined  using  codes  such  as  either  LOWTRAN 
or  MODTRAN.  Errors  in  the  codes  can  thus  lead  to  incorrect 
system  design.  Incorrect  designs  are  either  lacking  in 
sensitivity  or  have  excess  sensitivity  resulting  from 
overdesign.  Either  is  an  inefficient  use  of  resources. 

The  intent  of  this  thesis  is  to  determine  the  level  of  error 
in  these  codes  due  to  modeling.  This  characterization  can  be 
used  to  isolate  errors  in  both  sensor  design  and  data  analysis 
to  allow  for  more  efficient  utilization  of  scarce  remote 
sensing  resources. 
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2.  BACKGROUND 


Performance  of  this  thesis  required  an  understanding  of 
three  major  areas,  (A)  atmospheric  transmission  and  radiance 
modeling  programs  or  codes,  (B)  data  used  as  reference  against 
which  the  codes  will  be  compared,  and  (C)  the  statistical 
tools  and  terminology  used  in  comparing  the  code  output  with 
these  reference  data.  This  discussion  outlines  these  three 
major  areas,  their  underlying  theory,  and  the  procedures  used 
in  performing  the  work  associated  with  each  area. 

2.1.  Atmospheric  Transmission 

The  propagation  of  energy  is  central  to  remote  sensing. 
Remote  sensing  can  be  taken  to  imply  any  sensing  and  analysis 
of  events  while  not  in  direct  proximity  to  them.  This  thesis 
shall  employ  a  convention  that  defines  remote  sensing  in  a 
macroscopic  sense;  viewing  of  large  scale  events  over 
relatively  large  distances.  This  is  best  defined  by  using  the 
common  examples  of  airborne  or  satellite  sensing,  although 
ground  based  systems  can  also  be  included.  In  addition,  this 
thesis  focuses  on  remote  sensing  performed  in  the  part  of  the 
electromagnetic  spectrum  that  can  be  optically  imaged.  This 
spans  the  region  from  0.2  to  25  microns  (ultraviolet  through 
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deep  infrared) .  Energy  sensed  in  this  region  is  dominated  by 
two  sources,  propagation  of  the  energy  from  the  sun  into  the 
earth's  atmosphere  (primarily  in  the  visible  region)  and 
thermal  emissions  from  either  the  target  or  the  earth 
background  (dominant  in  the  long  wave  infrared  region) . 

Energy  passing  through  the  atmosphere  is  affected 
primarily  by  atmospheric  scattering  and  absorption  (Lillesand 
and  Kiefer,  1987)  \  Scattering  is  largely  dependent  on  the 
size  of  the  atmospheric  constituents  (molecules,  particles,  and 
aerosols)  compared  to  the  wavelength  of  the  energy.  Particles 
with  a  much  smaller  diameter  than  the  wavelength  of  the  energy 
cause  Rayleigh  scatter.  Rayleigh  scatter  is  inversely 
proportional  to  the  fourth  power  of  the  wavelength,  so  shorter 
wavelengths  are  scattered  much  more  than  longer  wavelengths. 
This  form  of  scattering  manifests  itself  in  the  blue  color  of 
the  sky.  Particles  with  diameters  approximately  the  same 
diameter  as  the  wavelength  result  in  Mie  scattering.  Mie 
scattering,  primarily  caused  by  dust  and  water  vapor,  effects 
longer  wavelengths  than  Rayleigh  scatter.  While  both  Rayleigh 
and  Mie  scatter  have  fairly  well  defined  scattering  patterns, 
not  all  scattering  is  so  well  behaved.  Nonselective  scatter, 
which  results  when  the  particle  diameter  is  larger  than  the 
wavelength  of  the  energy,  scatters  all  wavelengths  in  random 
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patterns. 

Absorption  is  the  other  primary  source  of  loss  in  energy 
transfer  through  the  atmosphere.  Absorption  is  caused  by  the 
interaction  of  energy  with  the  constituent  gases  of  the 
atmosphere;  including  water  vapor,  ozone,  and  carbon  dioxide. 
The  characteristic  energy  absorption  pattern  of  each  of  these 
gases  results  in  blocking  bands  (where  energy  is  absorbed)  and 
window  bands  (where  energy  passes  through  with  minimal 
absorption  effects) .  Absorption  is  an  inherently  wavelength- 
dependent  phenomenon. 

The  effects  of  atmospheric  transmission  losses  may  be 
seen  by  examining  Figure  1.  The  exoatmospheric  source 
function,  solar  radiation  reaching  the  earth's  outer 
atmosphere,  is  the  dotted  line.  The  irradiance  reaching  the 
sensor  is  shown  as  the  solid  line.  The  difference  between  the 
two  lines  is  caused  by  the  transmission  losses  described  in 
Figure  2.  Major  contributors  to  this  transmission  curve  are 
shown  in  Figures  3-6. 

Since  the  object  of  this  thesis  is  to  examine  solely  the 
effects  of  the  atmosphere  on  energy  propagation,  care  must  be 
taken  to  ensure  that  other  errors  are  not  induced.  Remote 
sensing  often  looks  at  the  reflectance  or  emission  of  a 
target.  This  thesis  avoids  the  uncertainties  of  that  problem 
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(eg.  specular  vs  lambertian  reflection,  wavelength  dependent 
emission,  effects  of  background)  by  looking  directly  at  the 
sun.  The  sun  may  be  approximated  as  an  exoatmospheric 
blackbody  with  a  temperature  of  ~6000°K.  Both  LOWTRAN  and 
MODTRAN  have  detailed  solar  source  models  which  will  make  this 
assumption  unnecessary.  If  this  is  taken  as  direct 
illumination  of  the  focal  plane,  the  only  other  sources  of 
energy  are  path  radiance  due  to  the  atmosphere  having  a  finite 
thermal  temperature  and  forward  scatter  of  solar  radiance. 

The  transmission  through  the  atmosphere  may  be  written 
using  Beer's  Law: 

Tp-e-P*  (1) 

where  Tg  is  the  transmission  through  the  path  length  z  and  R 
is  the  loss  coefficient  representing  losses  due  to  both 
scattering  and  absorption.  This  equation  is  solved  for  each 
homogeneous  path  of  length  z,  and  for  each  wavelength  of 
interest,  to  find  a  point  to  point  transmission  for  a  given 
wavelength.  A  homogeneous  path  is  defined  as  a  path  distance 
where  the  atmosphere  has  a  constant  temperature,  pressure,  and 
constituents . 

Radiance  at  a  sensor  then  may  be  approximately  written 
including  losses  as  (Schott,  1989)^: 


L-  (Egcosa  /n+ELp(l-x)  x^+EEgt^P  (0)  x" 


(2) 


where  L  is  the  radiance  received  at  the  sensor  (assuming  a 
standard  flat  focal  plane) .  Angles  are  defined  for  a 
hemisphere  (Figure  7) . 

This  term  may  be  more  easily  understood  by  viewing  it  as 
three  input  terms.  The  first  term  is  direct  solar  radiance: 

■^soJar  rad- (-^sCOSO  /u)  (3) 

where  is  the  exoatmospheric  solar  irradiance  outside  the 
atmosphere  on  a  plane  perpendicular  to  the  axis  of 
propagation,  a  is  the  angle  between  the  normal  to  the  surface 
and  the  sun,  and  is  the  optical  depth  through  the 

atmosphere  normal  to  the  earth  at  the  target. 

The  second  term  is  the  path  thermal  radiance: 

■^pac/j  zhernal  ^ 

where  Lj  is  the  radiance  associated  with  a  blackbody  at  the 
temperature  (T)  of  the  atmospheric  layer  and  (1-t)  is  the 
emissivity  of  the  layer,  and  t"  is  the  transmissivity  of  the 
path  from  the  layer  to  the  sun.  This  term  is  summed  over  the 
homogeneous  layers  from  the  point  of  emission  to  the  sensor. 
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The  third  term  is  the  forward  scattered  radiance: 

^forwazd  sent 

6(0)  is  the  angle  dependent  loss  coefficient  for  the  scattered 
emission,  and  and  r*'  are  the  transmissivity  of  the 
portions  of  the  path  viewed  prior  to  and  after  forward 
scatter.  This  term  is  summed  over  the  homogeneous  path 
distances  traveled  by  the  scattered  radiation. 

Thus,  the  energy  received  at  the  sensor  is  the  energy 
emitted  by  the  sun  corrected  for  the  angle  by  a  cosine  term, 
multiplied  by  the  transmissivity  of  the  atmosphere  (corrected 
for  the  angle) ,  and  the  thermal  path  radiance  summed  for  each 
layer  along  the  path  (corrected  for  the  emissivity  of  the 
path,  its  angle,  and  atmospheric  transmission) ,  and  forward 
scattered  energy  of  each  layer  along  the  path. 

This  equation,  like  equation  1,  must  be  solved  over 
homogeneous  paths,  z,  and  for  each  wavelength  of  interest. 
This  gives  rise  to  two  effects  that  must  be  corrected  for  in 
the  codes.  The  first  is  that  the  atmosphere  is  not  composed 
of  homogeneous  paths.  An  approximation  must  be  made  to 
generate  an  estimate  of  atmospheric  transmission,  either 
through  forced  establishment  of  homogeneous  layers  (LOKTRAN) , 
development  of  an  approximation  of  an  equivalent  homogenous 
layer  path  (LOKTRAIJ  5B  and  MODTRAN) ,  or  provision  to  allow  for 
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solution  of  equations  where  local  therinodynainic  equilibrium 
does  not  exist  (Fascode,  a  line  by  line  transmission  code  also 
developed  by  the  Air  Force  Geophysics  Laboratory,  (Smith, 
1978)^).  The  second  effect  is  known  as  curve  of  growth,  an 
effect  arising  from  the  pressure  broadening  of  absorption 
lines  as  a  function  of  optical  depth.  This  effect  is  not 
addressed  in  LOWTRAN  (since  it  is  a  20  cm’’  band  model  this 
effect  is  outside  its  effective  resolution  limits) ,  but  is 
addressed  in  LOWTRAN5B,  MODTRAN  and  Fascode. 

Neither  LOWTRAN  nor  MODTRAN  directly  solve  this  radiative 
transfer  equation.  Instead,  approximations  are  made  to  allow 
modeling  to  exist  within  the  framework  of  available 
information  and  analytic  tools. 

2.2.  LOWTRAN  Model 

The  LOWTRAN  model  is  a  low  resolution  (  20  cm*’  )  model 
originally  developed  by  the  Air  Force  Cambridge  Research 
Laboratory  (now  the  Air  Force  Geophysics  Laboratory)  in  1972 
(Selby  et  al,  1972)^.  Since  then,  considerable  improvements 
have  been  made,  expanding  the  code  to  include  radiance  in 
addition  to  transmission.  It  is  useful  to  review  the  form  and 
improvements  to  the  code  beginning  with  version  5.  LOWTRAN 
5(Knei2ys  et  al,  1980)^,  published  in  1980,  was  a  major 
software  revision  to  the  code  and  will  be  used  as  the  baseline 
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reference  for  the  rest  of  this  discussion,  as  the  extensions 
and  improvements (Kneizys  et  al,  1983)  (Kneizys  et  al,  1988)^'’’ 
incorporated  since  then  were  built  onto  version  5's  basic 
structure . 

LOWTRAN  5  calculates  atmospheric  transmission  based  upon 
molecular  band  absorption  and  scattering  (the  6  term  of 
Equation  1)  ,  aerosol  extinction,  and  molecular  continuum 
absorption.  Absorption  and  scattering  are  not  solved  directly 
for  the  Beer's  law  formulation  of  equation  1.  Instead,  a 
single  parameter  (absorption  coefficient.  S/d)  look-up  table 
is  employed  for  each  of  four  components.  These  components  are 
water  vapor,  ozone,  nitric  acid,  and  the  uniformly  mixed 
gasses  (COj,  NgO,  CH^,  CO,  Oj,  Ng) .  Each  band  is  represented  by 
a  single  parameter  generation  function  using  a  look-up 
parameter  based  upon  the  atmospheric  conditions  given  in  the 
inputs . 

Improvement  of  these  parameters  is  one  of  the  major 
changes  incorporated  in  the  most  recent  version,  LOWTRAN  7. 
The  need  to  upgrade  the  accuracy  in  the  parameters  chosen  for 
the  band  model  approach  was  outlined  in  Zachor  (1981)®.  This 
article  examined  the  accuracy  of  a  band  model  approach  to 
modeling  the  atmosphere.  Zachor 's  primary  conclusion  was  that 
the  single  parameter  solution  did  not,  and  could  not,  contain 
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enough  information  to  adequately  describe  the  atmosphere. 
This  was  particularly  true  where  multiple  gasses  were  to  be 
modeled  in  one  curve.  The  move  to  independent  curves  for  all 
of  the  major  constituent  gasses  in  LOWTRAN  7  addresses  these 
concerns. 

Using  a  method  described  by  Pierluissi  and  Tsai  (1987)’ 
the  uniformly  mixed  gases  are  now  treated  individually  rather 
than  as  a  group.  The  overall  formulation  is  similar  to  that 
of  the  earlier  LOWTRAN  versions  in  using  a  look-up  table 
rather  than  a  direct  solution  to  Equation  1.  The  effect  of 
incorporating  this  change  to  the  LOWTRAN  6  model  was  suggested 
and  reviewed  by  Cutten  (1986,  1988) Similarly,  the 
band  models  for  HjO  have  also  been  upgraded  (Pierluissi  et  al, 
1989)’^ 

The  continuum  model  for  water  vapor  was  upgraded  in 
LOWTRAN  6^  to  the  form  used  in  FASCODE  IB  (Clough  et  al, 
1981)’^.  This  provided  a  major  improvement  in  the  4. 5-5.0  lim 
region.  These  authors  note  that  there  are  still  great 
uncertainties  in  the  modeling  at  the  atmospheric  transmission 
windows  in  the  10  ^m  and  4  /xm  jcegions.  Improved  results  were 
achieved  in  the  4. 5-5.0  region. 

The  last  major  revision  to  the  physics  of  the  code  has 
been  in  the  treatment  of  scattering.  In  Version  5  all  energy 
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scattered  along  the  path  was  lost  and  no  energy  was  scattered 
back  in.  As  early  as  1980  Ben-Shalom  and  others  identified 
the  effects  that  this  would  have  over  long  path  distances^^. 
The  lack  of  conservative  scattering  is  estimated  to  cause 
errors  as  large  as  a  factor  of  two  under  some  viewing 
scenarios.  LOWTFIAN  6  introduced  a  single  scattering  model 
that  allowed  particles  to  scatter  once  before  being  lost. 
This  still  provided  results  significantly  different  from  those 
shown  using  fully  conservative  scattering  as  suggested  by  Ben- 
Shalom.  This  error  was  because  single  scattering  treated 
scattering  as  a  source  of  extinction  but  still  not  as  a  source 
of  radiance.  Fully  conservative  scattering  had  been 
implemented  in  Fascode  2^^,  but  this  proved  to  be  inaccurate 
as  well  since  now  all  energy  was  conserved,  leading  to  an 
overestimate  of  radiance. 

All  of  these  shortcomings  were  well  described  in  an 
article  by  Isaacs,  et  al  (1987)’^  in  which  they  describe  a 
multiple  scattering  modification  to  both  LOWTRAN  and  FASCODE 
which  greatly  reduces  errors  (claimed  to  be  less  than  20%  in 
all  viewing  geometries) .  Their  approach  is  to  model 
scattering  using  a  finite  stream  approach.  It  is  implemented 
in  LOWTRAN  7.  The  actual  results  depend  on  the  physical 
parameters  of  the  single  scattering  albedo  and  the  viewing 


20 


geometry  (representing  the  scattering  phase  function) . 

Another  modification  that  was  performed  to  the  LOVJTRAN 
model  was  the  inclusion  of  band  models  which  gave  true  5  cm-1 
accuracy.  (Standard  LOWTRAN  output  may  be  obtained  at  this 
resolution  but  it  is  interpolated  from  the  20  cm-1 
calculations.)  This  version,  known  as  LOWTRAN  5B,  developed 
by  Robertson  et  al  (1980)^^,  was  more  numerically  accurate 
than  the  interpolation  from  the  20  cm'^  output  of  LOWTRAN  since 
it  calculated  these  points.  It  also  predated  LOWTRAN  7's 
treatment  of  the  uniformly  mixed  gases  as  separate  parameters. 
However,  it  still  lacked  individual  band  models  for  all  of  the 
gasses  described  in  LOWTRAN  7  and  lack  of  multiple  scattering. 
Due  to  the  significant  differences  in  the  development  of  the 
band  model  parameters  it  will  be  discussed  in  detail  at  the 
beginning  of  the  MODTRAN  section. 

Both  LOWTRAN  and  MODTRAN  treat  the  atmosphere  as  layers 
above  the  earth  surface,  not  unlike  layers  of  an  onion  around 
its  core.  This  modeling  form  defines  both  how  energy  is 
propagated  through  the  atmosphere  (path  length)  and  how  the 
different  atmospheric  models  are  developed  with  their 
characteristic  components. 

The  LOWTRAN  5  manual  states  that  in  general  earth 
curvature  has  a  greater  effect  on  path  length  than  refraction, 


21 


except  at  angles  approaching  90°.  With  zenith  angles  close  to 
90°  the  increase  in  path  length  may  be  up  to  30%.  The  effect 
of  increasing  path  distances  is  to  increase  optical  depth. 
Optical  depth  refers  to  the  "thickness"  of  the  viewing 
atmosphere.  If  the  observer  is  looking  straight  up  the  depth 
is  one  atmosphere  as  each  layer  is  traversed  once  by  the 
energy  beam  along  the  shortest  (normal)  path.  As  the  viewing 
angle  is  changed  to  more  oblique  zenith  angles  each  layer  of 
the  atmosphere  has  a  longer  effective  depth  and  the  atmosphere 
appears  denser.  At  the  90  degrees  zenith  viewing  optical 
depth  is  approximately  39.65  atmospheres  (Bemporad  1907)’®. 

Refraction  within  the  atmosphere  is  governed  by  standard 
optical  principles,  principally  Snell's  law  for  the  transition 
of  energy  between  layers.  The  angle  of  refraction  from  one 
layer  to  another  is  defined  as: 

sin0-no  (Ro+Hj)  sin0(,/n(i?o+^)  (6) 


Where  G  is  the  zenith  angle,  n  is  the  index  of  refraction 
between  layers,  Rq  is  the  earth  radius,  H,  is  the  height  of  the 
lower  layer,  and  z  is  the  layer  depth.  The  subscript  0  on  the 
sine  term  indicates  arrival.  The  term  for  effective  path 
length  between  layers  Zj  and  Zj^,  is  derived  (Figure  8) 
where  B  is  the  angle  subtended  at  the  center  of  the  earth  as: 
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sinP  ^./sinO  j 


(7) 


Version  7  was  modified  to  account  for  surface  heights  of 
other  than  sea  level.  (Previous  versions  enforced  a  uniform 
earth  radius  with  all  surface  points  at  sea  level.)  LOWTRAN 
now  generates  a  modified  version  of  the  Beer's  law  formula; 
most  significant  at  altitudes  up  to  six  kilometers.  This 
modification  was  suggested  by  Novoseller  (1987)"’’  and  fully 
outlined  by  Shettle  (1989)^°.  The  effect  of  this 
modification  is  to  compress  the  atmospheric  layers  so  that 
areas  where  ground  level  is  above  sea  level  no  longer  suffer 
what  amounted  to  truncation  of  the  lower  atmospheric  layers. 
It  should  be  noted  that  this  does  not  change  the  molecular 
content  at  the  upper  altitudes,  rather  it  effectively 
compresses  the  effects  of  lower  altitudes  into  a  single  term 
so  that  the  extinction  is  properly  calculated. 

Molecular  content  at  an  altitude  is  defined  within  the 
scope  of  the  layer  concept.  Molecular  quantities  for 
important  species  are  defined  within  the  layers  at  predefined 
levels  for  the  different  aerosol  models  available  in  LOWTRAN 
and  MODTRAN.  The  different  models  attempt  ^o  provide  a 
reasonable  approximation  to  conditions  within  a  specified 
geographic  area  or  type  when  measurements  of  atmospheric 
conditions  are  not  available.  Most  changes  in  these 
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conditions  are  seen  in  the  lower  atmosphere.  Above  these 
lower  levels  (altitudes  >~10  km)  the  atmosphere  is  fairly 
uniform  globally  and  seasonally,  in  the  absence  of  volcanic 
activity. 

In  summary,  the  improvements  to  LOWTRAN  focus  on  the 
addition  of  multiple  scattering  and  separate  band  models  for 
the  most  Important  atmospheric  gasses.  This  is  not  to 
discount  other  major  improvements  such  as  the  revised  and 
highly  upgraded  extraterrestrial  source  function,  aerosol 
models,  and  cirrus  cloud  models;  it  simply  points  out  the 
major  changes  that  justify  a  review  of  LOWTRAN  capabilities 
now  that  more  accurate  results  can  be  expected. 

2.3.  MODTRAN 

MODTRAN  can  be  viewed  as  an  extension  of  the  band  model 
approach  used  in  LOWTRAN  5B  with  the  updated  aerosol  and 
scattering  functions  of  LOWTRAN  7.  Because  of  its  close 
similarities  to  MODTRAN,  LOWTRAN  5B  is  included  in  this 
section  and  will  be  discussed  first. 

LOWTRAN  5B  was  the  first  of  the  moderate  resolution 
codes.  It  retained  the  same  aerosol  models  as  LOWTRAN  5,  but 
utilized  a  significantly  different  method  of  developing  tlie 
band  model  parameters. 

It  included,  besides  the  increased  spectral  resolution 
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discussed  earlier,  addition  of  temperature  dependence  of  the 
molecular  absorption  coefficients,  and  use  of  a  multi 
parameter,  DctDpler-Lorentz  band  model. 

The  basic  forms  of  LOWTRAN  have  their  band  model 
coefficients  defined  for  an  absorption  coefficient,  S/d, 
defined  for  standard  temperature  and  pressure  (STP;  1  atm, 
273K) .  LOWTRAN  5B  calculations  for  the  coefficients  are  based 
on  multiple  temperatures  (200,  225,  250,  275,  and  300K)  and 
line  density,  1/d,  as  well  as  the  absorption  coefficient  S/d. 
Curve  of  growth  effects  are  included  through  the  combined 
Doppler-Lorentz  band  model.  The  Curtis-Godson  approximation 
is  used  for  multilayer  transmittance  calculations,  replacing 
an  inhomogeneous  path  with  a  homogeneous  one  by  using  average 
values  for  the  various  band  model  parameters. 

•'MODTRAN^^  is  a  moderate  resolution  model  and  computer 
code  used  to  predict  atmospheric  transmittance  and  background 
radiance  in  the  microwave,  infrared,  visible,  and  near 
ultraviolet  spectral  regions  (0  to  50,000  cm*’  or  0.2  micron 
to  infinity) .  The  code  maintains  complete  compatibility  with 
LOWTRAN  7,  specifically,  MODTRAN  retains  all  of  the 
capabilities  of  the  LOVJTRAN  7  model.  Both  codes  contain  the 
same  six  built-in  model  atmospheres;  spherical  refractive 
geometry,  aerosol  models,  clouds  (water  and  ice),  rain 
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attenuation  and  options  to  calculate  single  scattered 
solar/lunar  radiance,  solar/lunar  irradiance  and  multiply 
scattered  thermal  and  solar  radiance. 

"The  MODTRAN  code  improves  LOWTRAN's  spectral  resolution 
from  20  to  2  cm’’  full  width/half  maximum  (fwhm)  with  an  option 
to  vary  the  resolution  between  2  and  50  cm’’  (fwhm) .  The  band 
model  parameters  were  formulated  from  the  HITRAN  line  atlas 
for  twelve  atmospheric  gasses:  H2O,  CO2,  O3,  N2O,  CO,  CH^,  O2, 
NO,  SO2,  NO2,  NHj,  and  HNO3.  These  parameters  were  calculated 
for  1  cm’’  bins  from  0-17900  cm’’  at  five  temperatures  from  200 
to  300K,  all  stored  on  an  external  data  file  which  is  accessed 
by  the  program  (for  the  region  from  17,900  to  50,000  cm*’  the 
program  defaults  to  the  LOWTRAN  7  band  parameters.) 

"The  transmittance  is  calculated  with  an -equivalent-width 
formulation  which  accounts  for  the  finite  spectral  width  of 
each  interval  and  the  number  of  lines  contained  within  the  one 
cm’’  bin.  LOWTRAN  7  uses  a  one-parameter  band  model 
(absorption  coefficient.  S/d)  plus  molecular  density  scaling 
functions,  while  MODTRAN  relies  on  three  temperature  dependent 
parameters,  an  absorption  coefficient  (S/d) ,  a  line  density 
parameter  (1/d),  and  an  average  line  width.  The  absorption 
due  to  line  centers,  within  the  one  cm’’  bin,  is  modeled 
separately  from  the  absorption  due  to  line  tails  (from  regions 
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whose  line  centers  are  outside  the  specified  fwhm  interval) . 
The  absorption  due  to  lines  within  each  bin  is  calculated  by 
integrating  over  a  Voigt  line  shape.  The  line  tail  parameters 
consist  of  line  contributions  within  ±25  cm'^  (of  the  center 
of  the  fwhm  region)  .  The  line  tail  absorption  coefficient 
band  model  parameters  are  determined  by  integrating  the 
Lorentz  line  shape  over  this  interval  of  ±25  cm'^*"^^ 

The  improvements  to  the  band  model  parameters  and 
supporting  data  structures  (eg  solar  model,  scattering  model, 
and  aerosol  models)  as  well  as  the  introduction  of  the  true 
2cm'’  MODTRAN  model  suggests  that  computer  models  may  now  be 
expected  to  accurately  model  field  data. 

2.4  Field  Data 

Determining  the  transmission  and  radiance  accuracy  of  the 
two  codes  requires  some  form  of  standard  against  which  the 
output  is  judged.  In  this  research  the  object  is  to  isolate 
out  most  possible  sources  of  error.  The  simplest  means  of 
.doing  this  would  be  to  obtain  the  reference  measurements  in  a 
laboratory.  However,  laboratory  measurements  cannot  fully 
duplicate  the  effects  on  transmission  and  radiance  in  the 
atmosphere  over  long  path  lengths. 

However,  field  measurements  are  not  without  their  own 
shortfalls.  The  accuracy  of  the  reference  measurements  can  be 
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thought  of  as  a  function  of  the  accuracy  of  the  meteorological 
assumptions  of  the  code  inputs,  the  accuracy  of  the  reference 
source,  and  the  accuracy  of  the  instrument  taking  the 
measurement.  Chief  among  these  difficulties  is  obtaining 
meteorological  data  to  fully  describe  the  atmospheric  path 
through  which  the  measurement  is  taken. 

Ideally,  meteorological  data  would  be  available  over  the 
whole  path  of  interest.  This  is  unrealistic  for  long  paths. 
In  the  absence  of  meteorological  data  throughout  the  path 
model  inputs  will  be  used  instead. 

An  additional  difficulty  is  obtaining  a  good  reference 
signal.  In  this  thesis  the  sun  will  be  used  as  the  reference 
signal.  The  exoatmospheric  source  function  of  LOWTRAN  and 
MODTRAN  has  recently  been  upgraded  in  all  wavebands  using  data 
from  numerous  sources^^;  it  will  be  assumed  to  be  accurate 
(maximum  discontinuity  error  is  reported  as  three  percent  in 
the  midwave  infrared) ,  ignoring  the  day  to  day  variations 
induced  by  sun  flares  and  other  disturbances. 

The  accuracy  of  the  field  measurements  is  defined  by  the 
accuracy  of  the  instrument.  The  measurements  used  as 
reference  in  this  thesis  were  taken  by  the  Air  Force 
Geophysics  Laboratory  Flying  Infrared  Signatures  Technology 
Aircraft  (FISTA) .  The  particular  mission  reference  is  Mission 
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8202  flown  on  18  November  1981.  The  FI STA  program  is  a  suite 
of  instruments  flown  on  a  NKC--135  aircraft.  The  particular 
instrument  used  for  this  experiment  was  interferometer  105. 
Measurements  were  taken  of  the  setting  sun  over  the  Pacific 
Ocean  from  altitudes  of  39,000  ft  (11.89  km)  and  29,000  ft 
(8.84  km);  the  majority  of  the  measurements  were  obtained  at 
the  lower  altitude.  These  heights  are  sufficient  to  allow  use 
of  the  atmosphere  models  as  inputs  since  most  meteorological 
effects  occur  at  lower  altitudes.  Additionally,  the  s..y  was 
observed  to  be  cloud  free,  indicating  a  stable  tropospheric 
viewing  path. 

Interferometer  105  is  a  Michelson  type  interferometer. 
Its  output  is  a  Fourier  transform  of  the  actual  spectrum, 
which  is  inverse  transformed  to  arrive  at  the  spectral  data. 
The  theory  of  Fourier  transform  spectroscopy  is  covered  in 
numerous  articles  and  books  (Huppi,  et  al,  1981  and  Walker, 
1979)2526  will  not  be  discussed  in  this  thesis.  The 
FISTA  instrument  follows  conventional  design  and  analysis 
practices  for  this  type  instrument  (Draft  AFGL  document, 
1990)27. 

Of  greater  interest  is  the  sensitivity  and  resolution  of 
the  instrument  and  the  method  of  background  calibration  and 
removal.  The  parameters  of  concern  for  this  instrument  (are 
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stated  below  (Mills,  1989)^®: 

Wavelength:  2.0  -  5.4  41m 

Field  of  View:  1.0  degree  circular 

Spectral  Resolution:  0.95  cm’’ 

Detector:  InSb  (77K) 

Scan  Time:  1  second 

Sensitivity:  5E-T0  W/cm^  )Ltm. 

Solar  signatures  will  be  significantly  (~  factor  of  10^-10®) 
above  the  instrument  sensitivity  for  the  solar  viewing  cases, 
errors  should  be  confined  to  removal  of  background  effects. 
This  sensitivity  margin  will  not  exist  for  the  background 
signatures;  this  level  of  sensitivity  may  not  be  sufficient 
for  viewing  background  radiation. 

Target  irradiance  is  provided  as  the  calibrated  data 
product  from  FISTA  missions.  The  calibration  involves 
correction  for  off-axis  optics  effects.  This  calibration  is 
performed  prior  to  release  of  the  data.  Each  annular  region 
is  described  by  a  weighing  factor.  The  annular  zones  are  then 
used  to  weight  the  contribution  to  target  radiance  using  the 
form  of  the  following  equation: 

T-(T^-gBj/(l-g) 
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T  =  target  signal 
Tjj  =  center  radiance 
Bjj  =  outer  radiance 

g  =  fractional  value  assigned  to  outer  component  (Figure  9) . 
Since  FISTA  data  has  been  extensively  used  and  no  problems 
have  been  noted  with  its  calibration,  it  will  be  assumed  that 
the  FISTA  data  can  be  safely  used  as  a  standard  against  which 
to  evaluate  the  codes. 

2.5.  Statistical  Tools 

Traditionally,  differences  between  the  codes  and  a 
standard  have  been  defined  as  percent  error  over  the  waveband 
of  interest.  An  excellent  example  of  this  type  of  review  is 
available  in  a  report  by  Cundiff  (1986)^’.  In  many  aspects 
this  report  is  a  direct  precursor  to  my  thesis,  examining  the 
earlier  versions  of  these  codes.  The  percent  difference 
measure  is  adequate  only  over  discrete  wavebands  or  for 
broadband  integrated  radiance. 

A  more  statistically  rigorous  method  for  determining 
accuracy  of  the  code  generated  estimates  is  obtained  by 
comparing  the  code  output  directly  with  the  field  data,  point 
for  point.  This  comparison  allows  use  of  both  the  chi-squared 
and  Kolmogorov-Smirnov  tests.  The  Kolmogorov-Smirnov  test  is 
more  appropriate  for  testing  two  data  sets  and  is  conservative 
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for  count  data  (as  opposed  to  cumulative  percent  data) .  It 
is  a  nonparametric  test  designed  specifically  to  test  if  two 
independent  samples  come  from  identically  distributed 
populations  (Daniel,  1990)^°.  By  using  this  test,  as  well  as 
focusing  on  percent  differences  in  specific  wavebands,  the 
results  of  this  thesis  give  a  statistically  sound 
determination  of  the  errors  associated  with  the  two  codes. 

The  Kolmogorov-Smirnov  test  is  particularly  applicable 
since  the  statistical  tools  typically  used  (such  as  analysis 
of  variance)  require  an  assumption  of  normal  distributions,  an 
assumption  which  cannot  be  made  in  these  data.  General  linear 
model  techniques  are  inappropriate  for  the  determination  of 
significant  error  in  the  match  between  the  codes  and  data 
since  the  object  is  to  determine  the  degree  of  error  between 
points,  not  the  ability  of  a  polynomial  model  (equation)  to 
fit  the  data. 

2.6.  Previous  Efforts 

Characterization  of  the  errors  associated  with 
atmospheric  codes  is  not  a  new  or  unique  topic.  Examples 
include  many  of  the  references  already  provided  as  well  as 
others  of  a  more  general  nature  (Royer,  1988)^^.  Interest  in 
the  midwave  IR  bands  over  long  path  lengths  limits  the  number 
of  published  reports  considerably.  Examples  of  these  are 
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Cutten's  publications’'’®. 

Rather  more  indicative  of  the  efforts  performed  in  this 
area  are  internal  memos  and  reports  like  Cundiff^®.  This 
report  used  the  same  FISTA  data  as  this  thesis  but  earlier 
versions  of  the  codes.  Overall  accuracy  as  well  as  band 
accuracy  were  examined  as  percent  error  figures.  Extracts 
from  this  report  are  provided  as  Appendix  1.  The  percent 
differences  showed  maximum  error  in  the  bands  where  the  models 
were  the  weakest,  as  well  as  when  the  zenith  angle  was 
increased.  The  charts  representing  these  errors  are  presented 
in  pages  15  and  17  of  that  report. 

Summary  results  as  presented  from  page  37  of  that  report 

are: 

1.  Errors  increased  as  zenith  angle  increased, 

2.  Errors  were  largest  in  bands  where  the  least  modeling 
had  been  performed,  eg  the  2. 1-2.4  micron  CH,^  band, 

3.  Percent  error  measures  provide  adequate  diagnostics 
only  for  identifying  gross  trends, 

4.  Percent  error  is  inadequate  to  determining  overall 
quality  of  the  codes.  An  example  of  this  is  shown  in  the 
2. 1-2.4  urn  band  at  91.7°  zenith.  In  this  graph  there  are 
clearly  large  errors  on  the  order  of  100  percent,  but  the 
average  error  is  shown  as  near  50  percent. 
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Note  that  both  FASCODE  and  LOWTRAN  show  extreme  errors  at  near 
90°  zenith. 

The  errors  near  90°  zenith  may  not  be  a  failure  irt  the 
physical  modeling  but  rather  in  the  implementation.  Richter 
(1985)^^  demonstrates  that  the  dip  in  radiance  near  this 
angle  is  an  artifact  of  the  single  scattering  model.  It  is 
thus  reasonable  to  expect  significantly  improved  performance 
from  EOWTRAN  7  and  MODTRAN  near  90°.  This  is  due  to  the 
change  in  modeling  from  single  scatter  to  multiple  scattering 
effects  in  LOWTRAN  7  and  MODTRAN. 

Examination  of  individual  bands  is  shown  in  Cutten's  work 
as  well  as  in  a  follow  on  (Lisowski,  1988)^^to  Cundiff's  work 
also  provided  in  Appendix  1.  This  work  examined  the  COg  band 
in  an  attempt  to  chose  a  band  model  for  use  in  a  data 
reduction  process. 

The  discussion  shows  how  the  computer  models  have  evolved 
and  describes  the  field  data  and  statistical  tools  chosen  for 
evaluation  of  the  computer  models.  The  next  section  will 
describe  the  approach  used  for  this  evaluation. 

3 .  APPROACH 

The  goal  of  this  thesis,  a  statistical  estimate  of  error 
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between  a  computer  model  (LOWTRAN  or  MODTRAN)  and  a  baseline 
measurement  (the  AFGL  data  set) ^  required  generation  of  both 
a  calibrated  field  data  set  and  a  model  data  set.  It  was  then 
necessary  to  make  any  corrections  necessary  to  ensure  that  the 
two  data  sets  were  equivalent  before  making  the  error 
analysis.  This  section  is  logically  divided  into  a  discussion 
of  these  three  activities. 

3.1.  Field  Data  Generation 

The  data  utiliaed  in  this  thesis  was  obtained  on  18 
November  1982  by  the  FISTA  platform.  The  aircraft  flew  an 
orbit  off  of  the  Pacific  Coast;  with  March  AFB,  CA  the  takeoff 
and  landing  point.  During  the  orbit  the  setting  sun  was 
observed.  Matching  observations  of  the  background  sky  were 
obtained  on  each  orbit.  The  data  are  thus  in  pairs  with  the 
background  observations  oriented  180  degrees  from  the  sun 
(Figure  lOj.  The  physical  parameters  for  each  observation  pair 
are  provided  in  Table  1. 

The  AFGL  data  set  arrived  as  a  VAX  backup  tape. 
Unfortunately,  the  files  on  the  tape  were  in  a  format  unique 
to  the  AFGL  file  storage  system.  The  most  efficient  means  of 
reading  the  files  was  to  dump  the  files  in  both  word  and 
decimal  format.  These  files  were  then  edited  using  a  text 
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TABLE  1 

Air  Force  Geophysics  Laboratory  (AFGL) 

Flying  Infrared  Signatures  Technology  Aircraft  (FISTA) 


Mission  8202  18  November  1982 

Interferometer  105  (  1799.26-5199.15  wavenumbers) 

Run 

Pair 

GMT 

North 
Lat . 

West 

Long. 

Alt 

(km) 

Solar 

Elevation 

(Declin¬ 

ation) 

Solar 
Azimuth 
(“  East 
of 

North) 

1/2 

21.980 

34.67 

121.63 

11.887 

61.96 

214.29 

3/4 

22. 195 

34.96 

121.86 

11.887 

63.63 

216.94 

5/6 

22.428 

35.21 

121.96 

8.839 

65.55 

219.85 

7/8 

22.620 

34.96 

121.80 

8.839 

67.00 

222.49 

9/10 

22.859 

34.80 

121.60 

8.839 

69.05 

225.61 

11/12 

23.100 

34.83 

121.84 

8.839 

71.08 

228.21 

13/14 

23.286 

35.00 

12? .95 

8.839 

72.86 

230.17 

17/18 

23.736 

34.68 

121.59 

8.839 

77.31 

235.25 

19/20 

23.921 

34.86 

121.78 

8.839 

79.18 

236.94 

21/22 

0.069 

34.84 

121.63 

8.839 

80.81 

238.48 

23/24 

0.338 

35.00 

121.57 

8.839 

83.78 

241.03 

25/26 

0.530 

34.65 

121.81 

8.839 

85.53 

242.66 

27/28 

0.687 

34 .88 

121.84 

8.839 

87.34 

244.03 

29/30 

0.824 

34.82 

121.74 

8.839 

88.92 

245.29 

31/32 

0.976 

34.66 

121.84 

8.839 

90.48 

246.54 

33/34 

1.096 

34.68 

121.67 

8.839 

91.98 

247.65 

35/36 

1.205 

34.82 

121.68 

8.839 

93.26 

248.55 

37/36 

1.296 

34.76 

121.67 

8.839 

94.30 

249.31 
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Table  1  Notes: 

1.  60  dB  Neutral  Density  Filter  used  in  all  runs  except  37 

2.  Position  Locations  are  for  the  midpoint  of  each  data  run 

3.  Odd  numbered  runs  are  solar  observations,  even  numbered 
runs  are  background  observations. 

4.  Even  numbered  runs  use  position  information  generated 

during  matching  odd  run. 
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editor  to  confirm  the  text  information  in  the  file  header 
(already  available  in  text  format  in  the  information  provided 
by  AFGL  with  the  VAX  backup  tape,  Appendix  2) ,  and  edit  the 
data  into  a  numeric  format  suitable  for  efficient  processing 
on  the  computer.  A  sample  text  output  is  shown  below. 

Table  2 

Sample  Output  of  AFGL  Mission  8202  File  Header  Run  01 
Scale  Factor:  0.95367E-06 

Frequency  Spacing:  1.92847 

Beginning  Frequency:  1799.26 

Data  values  in  the  array: 1765 
AFGL  Mission  Designator:  MISSION  8202 
Archival  Tape  Designator : FILE  e:n00010  101.MIS202.  NUU 

FROM  R395-398/R491 
Notes:  SUN  ND  IN 

The  numeric  files  were  processed  using  the  Minitab 
statistical  spread  sheet  program.  The  routine  written  using 
this  software  read  in  the  decimal  numeric  file,  corrected  the 
order  of  the  data  (the  files  should  read  left  to  right  but  the 
numeric  file  was  inverse  order,  right  to  left),  stacked  the 
data  into  a  single  file  (the  AFGL  tape  contained  numeric  data 
in  two  records  which  together  formed  a  composite  numeric  data 
set) ,  and  then  ordered  the  data  in  column  format.  The  output 


of  this  process  was  a  file  with  the  raw  calibrated  counts, 
irradiance  in  watts /cin2 -cm,  the  wavenumber  for  the  particular 
irradiance  data  point  in  inverse  centimeters  and  the 
equivalent  wavelength  in  microns. 

Files  from  the  first  run  are  provided  as  an  example.  The 
decimal  numeric  file  is  Appendix  3.  The  resulting  data  file 
is  Appendix  4.  Due  to  its  length  (1768  data  points,  over  30 
pages  of  output)  only  the  first  and  last  pages  of  the  file  are 
provided  in  Appendix  4. 

The  relative  calibrated  irradiance  is  found  by 
multiplying  the  raw  count  number  (eg.  2026  for  the  first 
point)  by  the  scale  factor  (.95367E-06  for  this  file).  The 
effects  of  the  neutral  density  filter  utilized  in  this  mission 
for  solar  viewing  are  included  in  this  scale  factor.  The 
wavenumber  is  obtained  by  finding  the  offset  from  the  initial 
wavenumber.  The  initial  wavenumber  is  defined  as  1799.26  and 
the  step  interval  between  each  data  point  is  1.92847 
wavenumbers.  The  wavenumber  for  the  second  data  point  is  thus 
1799.26  +  1.92847  =  1801.19  wavenumbers. 

The  next  manipulation  was  to  derive  fully  calibrated 
data.  The  irradiance  values  generated  so  far  did  not  include 
any  correction  for  sensor  effects.  This  correction  was 
applied  by  multiplying  the  data  files  point-by-point  with 
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values  from  the  smoothed  correction  curve  (Figure  11) .  The 
smoothed  response  curve  is  derived  from  laboratory 
measurements  taken  in  a  Nitrogen  atmosphere.  The  effect  of 
this  curve  is  shown  in  the  plot  of  the  corrected  response 
curves  with  and  without  the  neutral  density  filter  (Figures  12 
and  13) .  Only  one  correction  curve  was  required.  (The 
neutral  density  filter  is  a  combination  of  a  metal  oxide 
filter  and  an  adjustment  to  the  instrument  electronics  and  has 
a  uniform  linear  effect  across  the  spectrum  of  interest. ^^) 
This  plot  closely  matches  other  laboratory  reference  curves 
for  this  and  other  AFGL  interferometers.^^  The  notch  visible 
in  both  response  curves  is  a  result  of  adding  a  correction 
factor  to  compensate  for  co^  gas  that  will  be  present  in  the 
aircraft  environment. 

It  was  necessary  to  determine  the  irradiance  in  units  of 
watts/cm^  -  micron.  A  simple  linear  interpolation  was  used  to 
convert  data  from  a  reference  to  inverse  centimeters  to  a 
reference  of  microns:  x  *  delta  wavelength/delta  wavenumber 
/  a  unit  correction  factor  =  y,  where  x  is  irradiance 
(wavenumbers)  and  y  is  irradiance  (microns) .  As  example,  for 
the  second  point,  .1.06E-06  *  (.0154/5)  /  lE-04  =  3.44E-04. 

One  final  manipulation  was  required  to  place  the  data  in 
the  format  desired  for  analysis.  The  AFGL  data  had  a  higher 
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spectral  resolution  than  LOWTRAN  could  output.  Data  to 
correspond  to  the  five  inverse  centimeter  LOWTRAN  output  was 
generated  by  interpolating  the  AFGL  data  to  that  resolution. 
This  was  accomplished  using  a  shape  conserving  cubic  spline 
interpolation  routine  available  in  the  IMSL^^  software 
package.  The  resulting  data  file  for  use  in  the  statistical 
comparisons  is  shown  in  Appendix  5.  (All  comparisons  were 
referenced  to  wavelength  values  in  microns.) 

3.2.  Computer  Model  Data  Generation 

Files  were  generated  using  both  LOWTRAN  7  and  MODTRAN  for 
direct  solar  illumination  and  radiance  with  multiple 
scattering  viewing  conditions.  The  physical  parameters  used 
were  those  recorded  by  AFGL  during  the  field  data  collection 
(Table  1) .  The  mission  data  sheet  did  not  provide  a  separate 
set  of  conditions  for  background  viewing.  Thus,  the  same 
physical  conditions  were  used  as  inputs  to  generate  both 
direct  solar  irradiation  files  and  the  paired  radiance  with 
multiple  scattering  files. 

An  annotation  on  the  flight  record  stated  that  no  clouds 
were  visible  in  the  viewing  path.  Due  to  the  extremely  long 
path  lengths  (ex.  183  km  for  Run  01  and  982  km  for  Run  29)  it 
was  necessary  to  utilize  the  standard  atmospheres  and  aerosol 
models  built  into  the  codes.  The  most  applicable  standard 
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atmosphere  was  the  tropical  atmosphere;.  All  of  the  other 
LOWTRAN  7  and  MODTRAN  model  atmospheres  were  tailored  to 
conditions  of  more  northerly  latitudes.  Table  2  shows  that 
similar  results  were  obtained  using  midlatitude  summer 
conditions  (the  next  closest  match  to  tropical  viewing 
conditions)  as  the  atmospheric  model. 

Several  different  aerosol  models  were  examined  as  inputs. 
The  most  logical  choice  was  the  Maritime  23  km  visibility. 
The  Navy  Maritime  model  was  not  used  since  it  required 
knowledge  of  surface  conditions  in  the  viewing  area  during  the 
preceding  24  hours.  This  data  was  not  available.  Analysis  of 
Table  3  indicates  that  the  most  critical  factor  regarding 
aerosols  in  the  computer  models  is  simply  their  inclusion  in 
the  model  inputs.  The  particular  choice  of  aerosol  model  is 
of  lesser  importance. 

The  data  for  the  LOWTRAN  files  were  generated  using  the 
IBM  Personal  Computer  (PC)  based  version  of  LOWTRAN.  Table  4 
shows  that  similar  results  were  obtained  for  both  the  PC  and 
mainframe  versions.  The  PC  data  was  used  for  comparison  as  it 
provided  the  most  rigorous  test  and  had  been  subject  to  the 
least  review  at  the  time  of  my  thesis. 

All  of  the  AFGL  data  was  in  units  of  irradiance  at  the 
sensor.  The  output  from  the  direct  solar  illumination  models 
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Table  3 

Effects  of  Aerosol  Models  Choice  on  Radiance  Values 
(Radiance  values  in  Watts/cm^-ster-cin) 


Tropical  Atmosphere 

Midlatitude 

Summer 

wave- 

number 

No 

Aerosols 

2  3  km 

Rural 

5  km 

Rural 

23  km 
Maritime 

23  km 

Maritime 

1800 

3.67E-09 

3 .75E-09 

3.75E-09 

3.64E-09 

3.75E-09 

2050 

7.77E-10 

8.46E-10 

8.46E-10 

1.02E-09 

8.46E-10 

2300 

4 . 78E-09 

4 .78E-09 

4.78E-09 

4.46E-09 

4.78E-09 

2550 

2.80E-11 

1 . 53E--10 

1. 53E-09 

1.57E-09 

1.53E-09 

2800 

1.51E-11 

1.77E-10 

1.77E-10 

1.79E-10 

1.77E-10 

3050 

1.33E-11 

2 . 14E-10 

2.14E-10 

2.14E-10 

2.14E-10 

3300 

1.24E-1I 

2.70E-10 

2.70E-10 

2.70E-10 

2.70E-10 

3550 

1.39E-11 

2.57E-10 

2.57E-10 

2.57E-10 

2.57E-10 

3800 

2.34E-11 

4 .37E-10 

4.37E-10 

4.37E-10 

4.37E-10 

4050 

4.06E-11 

7.11E-10 

7.11E-10 

7.11E-10 

7.I1E-10 

4300 

4.48E-11 

7.49E-10 

7.49E-10 

7.49E-10 

7.49E-10 

4550 

7.03E-11 

1.12E-09 

1.13E-09 

1.13E-09 

1.13E-09 

4800 

1.05E-10 

1.61E-09 

1.61E-09 

1.61E-09 

1.61E-09 

5050 

1.49E-10 

2.17E-09 

2.18E-09 

2.18E-09 

2.18E-09 
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Table  4 

Comparison  of  LOWTRAN  PC  and  Mainframe  Results 


Irradiance  (Watts/cm^-cm) 

wavenumber 

PC 

Mainframe 

%  Difference 

1800 

7.34E-07 

7.14E-07 

2.88 

2050 

9.64E-07 

9.37E-07 

2.88 

2300 

5.33E-09 

5.18E-09 

2.89 

2550 

1.55E-06 

1.51E-06 

2.65 

2800 

1.74E-06 

1.69E-06 

2 . 96 

3050 

1.91E-06 

1.85E-06 

3.24 

3300 

2.23E-06 

2.17E-06 

2.76 

3550 

1.74E-06 

1.69E-06 

2.96 

3800 

2.44E-06 

2.38E-06 

2.52 

4050 

3.22E-06 

3 .22E-06 

0.00 

4300 

2.88E-06 

2.79E-06 

3.23 

4550 

3.60E-06 

3.50E-06 

2.86 

4800 

4 . 36E-06 

4.24E-06 

2.83 

5050 

5.01E-06 

4 .87E-06 

2.87 

Mean  Percent  Difference  =  2.68 
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was  in  irradiance  units  but  the  output  of  the  radiance  with 
multiple  scattering  model  is  in  radiance  units.  Since  all  of 
the  AFGL  data  is  in  irradiance  units  it  was  necessary  to 
convert  the  output  of  the  radiance  models  to  irradiance  units. 
Using  the  formula  outlined  in  Slater (1980)^^  the  conversion 
from  radiance  to  irradiance  is  defined  by  knowing  the  field  of 
view  of  the  sensor.  The  sensor  has  a  nominal  one  degree  field 
of  view.  The  radiance  values  may  be  converted  to  irradiance 
values  for  this  sensor  by  multiplying  by  2 . 39E-04  steradians. 
(It  was  learned  after  completion  of  the  thesis  that  AFGL  uses 
an  empirically  derived  value  of  2.5E-04  value  for  converting: 
values  from  this  particular  instrument.) 

3.3.  Data  Set  Equivalence 

Even  though  both  sets  of  data  (field  and  computer 
generated)  had  now  been  established  as  calibrated  sets  with 
equivalent  viewing  conditions  and  units  a  final  correction  was- 
found  necessary.  While  attempting  to  compare  the  background 
files  it  was  found  that  two  sensor  effects  still  existed  that 
required  correction. 

The  first  effect  was  that  the  computer  models  were 
predicting  results  below  the  noise  floor  of  the  instrument. 
While  not  strictly  a  calibration  issue,  the  effect  was  the 
same.  Predictions  below  the  noise  floor  have  no  meaning  since 


51 


they  cannot  be  obtained  by  the  sensor.  The  correction  for 
this  issue  was  to  set  a  noise  floor  for  the  computer  values 
equal  to  the  stated  noise  floor  of  the  instrument.  Any 
predicted  values  below  this  level  were  adjusted  up  by  the 
value  of  the  noise  floor. 

This  adjustment  allowed  inclusion  of  points  that 
otherwise  would  have  had  to  have  been  discarded  from  the 
analysis.  While  it  does  not  allow  for  full  evaluation  of  the 
radiance  portion  of  the  models  it  does  allow  their  evaluation 
within  the  bounds  imposed  by  the  instrument. 

The  second  effect  was  apparent  blackbody  radiation 
affecting  data  in  the  longer  (5.5  micron)  end  of  the  spectrum. 
This  radiation  is  a  result  of  any  small  temperature  shift  in 
the  interferometer  sensor  system  between  the  times  of  the  cold 
calibration  and  the  field  collection.  This  effect  is  also 
believed  to  include  emission  from  the  telescope  barrel  and  the 
instrument  optics.  This  effect  that  is  not  removed  using  the 
smoothed  correction  curve^®.  (The  smoothed  correction  curve 
is  developed  from  the  cold  calibration  data.)  The  effect  was 
on  the  order  of  7E-09  Watts/cm^-micron  at  its  maximum  (near 
the  5.5  micron  region). 

The  correction  was  to  add  greybody  radiation  as  if  the 
optics  had  been  a  blackbody  source.  A  series  of  curves  was 
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plotted  against  the  data-  These  curves  utilized  blackbody 
temperatures  from  270  to  285K  and  emissivities  from  .1  to  .2. 

A  curve  from  a  blackbody  with  emissivity  of  15%  and 
temperature  of  285  degrees  Kelvin  was  found  to  provide  a  good 
fit  to  the  data  (Figures  14-17).  These  values  were  felt  to  be 
reasonable  since  the  optics  material,  CaFj  has  a  stated 
emissivity  of  10%  for  pure  material^’  and  the  temperature  is 
also  a  reasonable  approximation  for  a  transition  between  the 
outside  and  inside  of  the  aircraft  platform.  Figures  18  and 
19  provide  a  graphic  reference  for  the  similarities  between 
the  calibrated  MODTRAN  and  LOWTRAN  and  the  AFGL  data. 

4 .  RESULTS 

Two  sets  of  tests  were  made  on  the  data  sets.  The 
results  of  the  Kolmogorov-Smirnov  test  determine  if  the  field 
data  and  computer  data  follow  the  same  distribution.  The 
results  of  the  percent  difference  test  provide  details  at 
specific  wavelengths  and  insight  into  error  trends. 

Equivalent  inputs  were  used  in  both  sets  of  tests.  For 
background  viewing  conditions  (even  numbered  AFGL  runs)  the 
comparison  took  the  form: 


53 


CO 

00 

CO 

00 

O) 

c> 

<i> 

6 

0 

d> 

U) 

o 

lO 

o 

o 

cvi 

csi 

ui 

(uoj3!iu-3uio/s))eM}  aoueipejjj 

54 


I 


with  285K  BB  and  .15  Emitter  - -  afgl"'® 


r  «> 


I  I"  I  J  I  •"  I  '1  'I  !■  t—  t  '  I - 1  I  I'  ■  I  r~r  I  'I  I  V  »'  I  *"■ 


00 

op 

00 

00 

o 

o 

o 

6 

o 

o 

ii 

+ 

o 

o 

o 

o 

o 

o 

U) 

o 

iq 

o 

o 

o 

o 

oi 

cj 

U) 

d 

(uoj3!Ui'3UJD/$ueM)  eoueipejji 

57 


microns 


MODTRAN  and  AFGL  Data 


microns 


(uojoiui-zuio/s)}eM)  aoueipejji 
59 


microns 


“  ^AFGL  ^RadiancewithMulzipleScactering  “ 

For  solar  viewing  conditions  (odd  numbered  AFGL  runs)  the 
comparison  took  the  form: 

“  ^AFGL  ^Dii  Solar  IIT  ^Rad  with  Mult  Scat  “  (IQ) 

(AFGL  solar  viewing  captures  at  the  instrument  both  direct 
solar  irradiance  and  irradiance  from  the  sky  background  with 
multiple  scattering  of  both  sky  and  solar  radiance) .  Thus, 
for  all  comparisons,  E„  is  the  measured  irradiance  (field 
data)  and  Ej.  is  the  equivalent  calculated  irradiance  (computer 
data)  at  the  instrument. 

4.1.  Kolmogorov-Smirnov  Analysis 

The  Kolmogorov-Smirnov  test  is  a  hypothesis  test  that  two 
sets  of  data  have  the  same  probability  density  function.  The 
form  of  the  hypothesis  is  a*^  follows  for  the  two  sided  case: 

Hq{x)  :  Fyix)  ~  F2ix)  for  all  X  (11) 

F^ix)  *F2{x)  for  at  least  one  X  (12) 

The  test  results  (Table  5)  are  generated  in  terms  of  the  two 
sided  test  statistic,  D,  and  are  shown  in  terms  of  the 
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Table  5 

Results  of  Kolmogorov-Smirnov  Test 


AFGL  Run  Number 

Probability  of  Greater  D,  Two-Sided  Test 

(Solar  Elevation) 

LOWTRAN  Comparison 

MODTRAN  Comparison 

1  (61.96°) 

0.0001 

0.0003 

2  (61.96°) 

0.0000 

0.0000 

3  (63.63°) 

0.0003 

0.0010 

4  (63.63°) 

0.0000 

0.0000 

5  (65.55°) 

0.0073 

0.0206 

6  (65.55°) 

0.0000 

0.0000 

7  (67.00°) 

0.0005 

0.0034 

8  (67.00°) 

0.0000 

0.0000 

9  (69.05°) 

0.0007 

0.0042 

10  (69.05°) 

0.0000 

0.0000 

11  (71.08°) 

0.0004 

0.0034 

12  (71.08°) 

0.0000 

0.0000 

13  (72.86°) 

0.0000 

0.0000 

14  (72.86°) 

0.0000 

0.0000 

17  (77.31°) 

0.1735 

0.1195 

18  (77.31°) 

0.0000 

0.0000 

19  (79.18°) 

0.0388 

0.0284 

20  (79.18°) 

0.0000 

0.0000 

21  (80.81°) 

0.3368 

0.1537 

22  (80.81°) 

0.0000 

0.0000 

23  (83.78°) 

0.0332 

0.0332 

24  (83.78°) 

0.0000 

0.0000 

25  (85.53°) 

0.0801 

0.0451 
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26  (85.53°) 

0 . 0000 

0.0000 

27  (87.34°) 

0.0206 

0.0061 

28  (87.34°) 

0.0000 

0.0000 

29  (88.92°) 

0.0000 

0.0000 

30  (88.92°) 

0.0000 

0.0000 

31  (90.48°) 

0.0000 

0.0000 

32  (90.48°) 

0.0000 

0.0000 

33  (91.98°) 

0.0000 

0.0000 

34  (91.98°) 

0.0000 

0.0000 

35  (93.26°) 

0.0000 

0.0000 

36  (93.26°) 

0.0000 

0.0000 

37  (94.30°) 

0.0000 

o.oopo 

Table  4  Notes: 

1.  Odd  Runs  are  Solar  Viewing,  Even  Runs  Background  Viewing. 
2i  Runs  1-4  are  at  11. 887kin,  all  others  are  at  8.839kTn. 
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resulting  probability.  Probability  is  defined  as  l~alpha 
where  alpha  is  the  risk  of  rejecting  the  null  hypothesis  when 
it  is  true.  Table  4  shows  that  for  most  of  the  runs  Hq  is 
true  with  a  probability  greater  than  97%.  It  is  thus 
concluded  that  the  two  sets  of  data,  field  and  computer,  have 
the  same  probability  distribution  function.  From  this  it  can 
be  inferred  that  they  have  statistically  the  same  shape  and 
that  the  computer  models  provide  a  good  representation  of  the 
observed  irradiances. 

4.2.  Percent  Difference  Analysis 

The  percent  difference  was  defined  as: 


A%(X) 


(13) 


where  Ep(lambda)  is  the  calculated  value  at  a  given  wavelength 
and  E^( lambda)  is  the  measured  value  at  a  given  wavelength. 
A  positive  percent  difference  represents  an  under  estimate  of 
the  irradiance  by  the  models. 

Singularly  large  differences  at  a  given  wavelength  can  be 
generated  using  equation  13,  especially  in  the  blocking  bands 
of  HjO  at  2.7  microns  and  CO^  at  4.3  microns.  These  large 
differences  arise  from  small  differences  between  the  actual 
concentrations  of  these  atmospheric  gasses  and  the 
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concentrations  predicted  by  the  computer  codes.  In  the 
blocking  bands  small;  changes  in  concentration  result  in 
relatively  large  changes  in  transmissivity.  These  large 
changes  in  transmissivity  of  the  atmosphere  have  large  effects 
on  the  irradiance  at  the  sensor.  Thus,  it  is  possible  to 
arrive  at  large  percent  changes  in  irradiance  from  errors  in 
the  calculated  data  that  result  not  from  errors  in  the  model 
itself,  but  from  the  atmospheric  parameters  input  to  the 
model . 

Therefore,  a  limit  was  applied  to  the  results  of  the 
percent  difference  calculations.  This  limit  sets  results  from 
equation  13  that  exceed  a  percent  difference  of  ±200%  to 
±200%.  This  was  done  so  that  singularly  large  differences  at 
a  given  wavelength  would  not  impede  generation  of  useful  plots 
or  analysis  at  other  wavelengths.  The  number  of  "limited" 
points  can  also  be  used  as  a  diagnostic  tool  to  evaluate  the 
level  of  error  in  a  given  comparison;  runs  with  larger  numbers 
of  "limited"  points  have  greater  overall  error  (Table  6)  .  A 
tabular  summary  of  the  statistics  is  provided  (Table  7). 
Plots  of  these  values  are  provided  as  Figures  20-29. 

The  choice  to  "limit"  the  data  was  made  over  the 
alternative  of  declaring  data  in  these  regions  invalid,  as  in 
Cundiff .  Either  approach  has  merit.  Since  the  intent  of  this 
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Table  6 

Summary  Results  from  "limiting" 
Percent  Difference  Calculations 


#  of  LOWTRAN 
points 
"limited" 


%  of  LOWTRAN 
points 
"limited" 

#  Of  MODTRAN 
points 
"limited" 

%  of  MODTRAN 
points 
"limited" 

2.35 

14 

2.06 

15.63 

105 

15.49 

2.51 

16 

2.35 

17.11 

116 

17.11 

3.24 

30 

4.43 

11.50 

80 

11.80 

3.24 

27 

3.98 

10.91 

73 

10.76 

3.39 

32 

4.72 

11.21 

74 

10.91 

3.98 

37 

5.46 

8.41 

59 

8.70 

3.98 

36 

5.31 

7.52 

53 

7.82 

6.05 

43 

6.34 

25.66 

172 

25.37 

5.75 

44 

6.49 

4.57 

30 

4.42 

6.19 

55 

8.11 

3.98 

30 

4.42 

7.96 

62 

9.14 

1.77 

134 

19.76 

65 


25 

70 

10.32 

57 

11.06 

26 

12 

1.77 

127 

18.73 

27 

87 

12.83 

92 

13.57 

28 

6 

0.09 

123 

18.14 

29 

112 

16.52 

116 

17.11 

30 

9 

1.33 

59 

8.7 

31 

162 

23.89 

185 

27.29 

32 

7 

1.03 

57 

8.41 

33 

472 

69.62 

509 

75.07 

34 

13 

1.92 

49 

7.22 

35 

676 

99.71 

676 

99.71 

36 

25 

3.69 

71 

10.47 

37 

4 

0.59 

26 

3.83 

66 


Table  7 

Statistics  for  "limited"  Percent  Difference  Comparisons 


Statistics  of  "limited" 
LOWTRAN  Comparisons 


Statistics  of  "limited" 
MODTRAN  Comparisons 


AFGL 
Run  # 

Mean 

Median 

Std  Dev 

Mean 

Median 

Std  Dev 

1 

4.10 

5.90 

35.78 

3.87 

5.79 

33.35 

2 

-43.45 

-12.08 

91.13 

-43.66 

-11.66 

91.25 

3 

2.49 

4.86 

36.58 

1.67 

4.78 

34.96 

4 

-35.14 

6.88 

96.93 

-35.24 

7.17 

97.02 

5 

-2.07 

3.03 

44.14 

-5.43 

3.34 

47.61 

6 

-29.98 

-1.65 

86.28 

-30.27 

-2.01 

86.30 

7 

-3.10 

3.22 

45.11 

39.71 

51.52 

52.95 

8 

-31.96 

-4.21 

84.95 

-31.66 

-4.69 

84.97 

9 

-2.72 

3.27 

45.46 

-6.72 

3.06 

50.04 

10 

-31,87 

-3.77 

88.53 

-32.25 

-4.40 

88.67 

11 

-3.97 

3.17 

46.72 

-7.96 

3.70 

51.64 

12 

-20.46 

4.48 

82.40 

-20.86 

4.49 

82.67 

13 

0.12 

9.00 

49.15 

-2.25 

9.07 

52.95 

14 

-18.88 

-1.96 

81.79 

-19.34 

-2.41 

81.94 

17 

-18.21 

-7.49 

52.49 

-20.38 

-7.88 

53.11 

18 

-54.40 

-29.78 

107.96 

-55.05 

-30.87 

108.09 

19 

-22.65 

-11.79 

52.19 

-25.38 

-11.74 

53.39 

20 

-13.22 

-1.73 

73.84 

-13.94 

-2.04 

74.21 

21 

-15.24 

-3.58 

55.58 

-22.45 

-5.13 

58.42 

22 

-2.20 

16.74 

12.46 

-2.97 

15.55 

73.18 

23 

-12.83 

3.05 

62.14 

-15.53 

2.30 

63.38 

24 

5.44 

9.82 

63.81 

-41.35 

-6.26 

96.17 

67 


25 

-23.65 

-4.21 

67.43 

-27.44 

-5.70 

69 . 11 

26 

5.44 

6.19 

62.76 

-50.03 

-17.01 

94.06 

27 

-25.51 

0.38 

73.58 

-28.11 

-1.20 

73.43 

28 

1.96 

0.03 

58.41 

-62.24  : 

-39.11 

88.38 

29 

-27.10 

5.22 

82.59 

-31.61 

2.52 

81.97 

30 

27.81 

30.43 

55.75 

-22.78 

5.18 

82.63 

31 

-52.44 

-8.09 

91.59 

-65.19 

-18.61 

91.27 

32 

25.02 

31.67 

58.29 

-38.58 

84.97 

33 

-167.25 

-200.00 

57.77 

-169.68 

-200.00 

58 . 02 

34 

40.50 

54.28 

65.43 

2.47 

26.17 

82.43 

35 

-199.12 

-200.00 

16.28 

-199.12 

-200.00 

16.28 

36 

41.22 

70.84 

75.10 

3.10 

30.72 

90.13 

37 

66.82 

97.22 

50.25 

46.35 

71.09 

73.17 
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Number  of  MODTRAN 
"limited”  Points 
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AFGL  Run  Number 


Figure  23 
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LOWTRAN  Mean  Percent  Difference 
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thesis  was  to  evaluate  performance  over  the  full  ?  0-5.5 
micron  region,  "limiting”  the  data  provided  better  ins. 3ht 
into  the  errors  of  the  codes.  In  addition,  it  is  difficult  to 
decide  on  an  single  set  of  intervals  that  should  be  excluded 
as  invalid  due  to  curve  of  growth  effects. 

The  optimum  approach,  which  would  have  avoided  the 
uncertainty  in  band  exclusion  intervals  would  have  been  to 
perform  a  signal  to  noise  analysis.  Data  with  a  signal  to 
noise  ratio  of  less  than  three  could  be  regarded  as  having 
insufficient  strength  to  justify  analysis.  Since  the  AFGL 
data  was  provided  as  a  processed  signal  this  analysis  could 
not  be  performed  using  this  data  set. 

Plots  for  every  fourth  run  are  provided  as  Figures  30-65 
in  Appendix  8.  Analysis  of  these  plots  and  Table  7  provided 
the  following  observations-: 

4.2.1.  Median  errors  are  relatively  small  in  the  lower  zenith 
angles.  Errors  tend  to  stay  below  an  average  of  10%  at  zenith 
angles  less  than  90  degrees  (Runs  1-30),  as  example  Figure  30. 
Errors  are  well  above  10%  at  zenith  angles  equal  to  or  larger 
than  90  degrees  (Runs  31-37)  (Figure  37) .  The  results  from 
the  direct  solar  viewing  cases  indicate  that  the  models  for 
transmission  and  solar  irradiance  are  essentially  correct 
(Figures  66  and  68) ,  except  for  the  regions  noted  in 
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paragraphs  4.2.5  and  4.2.6.  The  plots  of  median  error  for 
background  viewing  indicate  a  slightly  larger  error  band  at 
low  zenith  angles,  but  without  the  extreme  excursions  at  high 
zenith  angles  (Figures  67  and  69) .  However,  this  result  may 
be  due  to  the  large  number  of  points  affected  by  the 
correction  for  the  sensor  noise  floor  and  greybody  optics. 
4.2.2.  Errors  occur  largely  in  the  regions  of  the  primary 
atmospheric  absorbers. 


Wavelength  Region 
(microns) 

2.0 

2.75 

3.3 


Primary  Atmospheric 
Absorber 


H2O 


4.3 


4.8 


5.5 


HjO 


CH, 


CO, 


CO2  /  Oj 


H2O 


4.2.3.  The  large  percent  differences  in  the  major  blocking 
regions,  2.7  micron  H2O  and  4.3  micron  CO2,  are  numerical 
artifacts.  The  low  transmission  of  these  regions  (Figures  3 
&  4)  mean  that  relatively  small  changes  in  transmission 
(irradiance  at  the  sensor)  result  in  disproportionately  large 
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Figure  66 

MODTRAN  Odd  Run  Percent  DIffernce 
Plotted  Against  Degrees 
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LOWTRAN  Odd  Run  Percent  Difference 
Plotted  Against  Degrees 
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percent  difference  measurements. 

4.2.4.  The  increased  viewing  of  higher  density  regions  of  the 
atmosphere  results  in  broadened  absorption  lines.  Curve  of 
growth  effects  are  clearly  visible  in  the  results  of  the  2,15, 
3.3,  and  4.3  micron  regions  for  both  MODTRAN  and  LOWTRAN. 
This  causes  the  bands  of  large  percent  difference  in  these 
regions  to  "grow"  in  width  as  the  zenith  angle  increases. 

4.2.5.  The  error  trend  in  the  2.0  -  2.4  micron  region  of  the 
solar  viewing  observations  appears  to  have  a  consistent  shape. 
This  shape  arrears  independent  of  zenith  angle  and  is  the 
similar  for  both  MODTRAN  and  LOWTRAN.  The  error  indicates 
that  the  models  are  predicting  more  solar  irradiance  than  is 
received  by  the  sensor  (Figures  18  and  19) .  Since  this  error 
does  not  change  with  angle  it  is  not  likely  to  be  a  result  of 
error  in  transmission,  which  will  change  with  view  angle.  The 
two  other  major  sources  of  error  are  errors  in  the  solar  model 
or  in  the  instrument  calibration.  Since  the  solar  model  was 
compiled  based  on  multiple  observations  taken  by  many 
observers  over  time  it  is  not  as  likely  to  contain  the  -20% 
error  indicated  in  the  results.  This  suggests  that  the 
calibration  curve  for  the  AFGL  instrument  is  in  error  in  this 
wavelength  region. 

4.2.6.  The  error  in  the  5.2  -  5.5  micron  region  is  angle 
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dependent.  This  indicates  a  possible  failing  in  the 
transmission  modeling.  The  primary  constituent  of  this  region 
is  H2O.  As  noted  in  the  discussion  (Section  2) ,  the  modeling 
of  the  continuum  transmission  of  HjO  continues  to  be  an  area 
of  research.  Errors  in  the  transmission  could  come  from  the 
water  content  provided  by  the  standard  model  and  from  the 
transmission  curve  for  water  vapor  continuum.  Since  errors 
exist  in  the  primary  water  band  at  2.7  microns  as  well  as  in 
this  band  it  is  suspected  that  the  error  in  the  5.2  -  5.5 
region  is  a  result  of  both  sources  of  error. 

4.2.7.  Standard  deviation  measurements  (Figures  28  &  29)  are 
a  poor  indicator  of  quality.  Normally  a  smaller  standard 
deviation  is  an  indication  of  better  data  since  it  means  there 
is  less  variance  in  the  data.  In  the  case  of  the  percent 
difference  measurements  this  is  a  false  indication  since  the 
reduction  of  variation  is  a  result  of  greater  "limiting”  of 
the  results.  The  runs  with  lowest  standard  deviation  are  also 
the  runs  with  the  greatest  number  of  "limited"  points.  These 
runs  are  at  the  highest  zenith  angles,  already  noted  as  the 
runs  with  the  largest  errors. 

4.2.8.  The  difference  in  results  in  LOWTRAN  and  MODTRAN 
results  in  the  low  zenith  angle  runs  where  errors  are 
essentially  zero  is  an  indication  that  the  2.68%  upward  bias 
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of  the  LOWTRAN  PC  generated  values  (shown  in  Table  4)  is  true 
bias  and  not  just  an  artifact.  This  would  indicate  an  upward 
bias  of  the  PC  based  LOWTRAN  results  compared  to  field  data. 
These  results  suggest  that  if  accuracy  is  critical  then  a 
mainframe  computer  should  be  used  to  perfora  the  calculations. 
(As  part  of  the  installation  of  MODTRAN  a  limited  comparison 
was  made  between  MODTRAN  and  LOWTRAN  7  results  when  both  codes 
were  run  on  the  same  mainframe.  This  limited  comparison 
showed  the  two  codes  generating,  results  that  agreed  much  more 
closely  than  the  results  shown  in  Table  4.  This  suggests  that 
the  bias  shown  in  Table  4  is  due  largely  to  the  computer,  hot 
the  model.) 

4.2.9.  Close  examination  of  Table  7  shown  the  following 
discrepancies  from  the  general  trends:  noted  above 

4;. 2. 9. a.  There-  is  a  difference  between  the  MODTRAN  and 
LOWTRAN  results  for  Runs  7i,  26,  28,  and  30.  The 
difference  for  Run  7  lacks  an  explanatjon.  The 
difference  for  runs  26,  28,  and  30  (Figures  45  and  (53  for 
Run  28)  may  he  due  to  the  increase  in  observed 
Irradiance.  For  runs  26  and  28  LOWTRAN  sees  a  lower 
value  due  to  its  coarser  structure.  It  thus  has  a 
greater  number  of  matches  enforced'  by  the  fit  to  the 
sensor  noise  floor.  By  run  30,  MODTRAN  has  sufficient 


irradiance  to  achieve  a  match  while  LOWTRi^  is  now  seeing 
just  sufficient  irradiance  to  be  above  the  noise  floor 
but  not  enough  to  achieve  a  match. 

4. 2.?.b.  Increasing  errors  for  both  MODTRAN  and  LOWTRAN 
are  observed  in  both  Runs  18  and  22.  This  is  due  to  an 
insufficient  match  to  the  data  by  the  greybody  and  sensor 
noise  corrections. 

4. 2.9.C.  Runs  33  and  35  exhibit  extreme  error  for  both 
LOWTRAN  and  MQDTRAN. 

5.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  results  of  my  thesis  strongly  parallel  the  results 
reported  by  cundiff?®  in  his  analysis.  As  with  Cundiff ,  errors 
increased'  as  zenith  angle  increased.  MODTRAN  results  closely 
mirror  his  results  for  LOWTR^  5B,  an  expected  result  since 
the  band  models  share  a  similar  formulation.  LOWTRAN  7  shows 
a  significant  improvement  over  LOWTRAN  6,  as  expected,  both 
models  show  improved  performance  in  the  4.3  fin  region, 
compared  to  the  results  obtained  by  Cundiff. 

Several  overall  conclusions  can  be  drawn  from  the 
analysis  and  results. 

5.1.  LOWTRAN  and  MODTRAN  give  similar  results.  Howeven,  the 
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bias  in  the  PG  generated  LOWTRAN  data  (most  clearly  indicated 
in  the  solar  viewing  cases  at  low  zenith  angles)  indicates 
that  the  mainframe  versions  of  the  codes  should  be  used 
whenever  practical  to  obtain  the  most  accurate  resultSi  In 
the  region  reviewed  in  this  thesis  no  preference  is  indicated 
for  either  model.  This  conclusion  may  be  due  to  the  smoothing 
effect  of  the  i^nterpolation  of  the  field  data  to  5  cm*’ 
resolution,  the  similar  smoothing  effect  obtained  in  the 
interpolation  of  LOWTiy^  20  cm*’  resolution  data  to  5  cm*’ 
data,  and  limiting  output  resolution  of  the  MODTRAN  model  to 
the  5  cm*’  resolution  of  LOWTRi^. 

5.2.  The  computer  models  have  error  results  within  10%  for 
moot  viewing  cases  at  zenith  angles  less  than  90  degrees 
(Figures  67  and  68) .  This  is  especially  true  where 
transmission  is  the  major  factor. 

5.3.  ‘Results  for  radiance  with  multiple  scattering  require 
further  investigation.  The  adjustments  for  sensor  noise  flOor 
and  greybody  radiation  precluded  full  analysis  of  this  part  of 
the  computer  models.  The  large  number  of  points  affected  by 
these  corrections  and  the  relative  magnitude  of  the  correction 
compared  to  the  predicted  values  leaves  doubt  whether  any 
observed  differences  are  derived  from  the  corrections  imposed 
for  analysis  Or  a  result  of  differences  between  the  models  and 
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real  world  conditions. 

5.4;.  The  inability  to  aodel  the  fine  structure  of  the  CH^ 
lines  in  the  3.3  micron  region  indicates  a  fundamental 
limitation  of  any  band  model  code.  If  specific  structural 
features  are  important  in  either  data  analysis  or  sensor 
design  then  a  line-by-line  code  such  as  Fascode  is  required. 
5.5.  Improvements  are  required  in  the  continuum  modeling  of 
HjO  in  the  2:.  1  and  5 .5  micron  regions .  These  improvements 
should  correct  the  errors  in  the  transmission  curves  in  these 
regions.  Since  the  errors  in  these  regions  could  also  be  a 
result  of  errors  in  the  atmospheric  input  parameters  to  the 
codes  it  is  recommended  that  research  into  this  area  utilize 
data  from  a  set  of  measurements  where  the  atmospheric  viewing 
path  is  closely  defined  both  in  geometry  (little  dr  no 
refraction)  and  in  atmospheric  constituents. 

Further  cdhclusions  require  acquisition  of  additional 
data  to  verify  the  radiance  with  multiple  scattering  models. 
The  sensor  used  to  obtain  this  data  must  have  sufficient 
sensitivity  to  conf  irm  radiance  modeling  results  in  the  region 
2.7  -  4.5  microns.  It  was  this  region  that  was  limited  by  the 
FISTA  instrument  noise  floor.  Calibration  shduld  include 
matching  of  the  response  in  the  HjO  regions  and  in  the  CO2 
region  at  4 . 3  microns  to  allow  discrimination  between  sensor 
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calibration  and  computer  model  errors. 

It  is  suggested  that  a  separate  set  of  measurements 
(similar  to  Cutten/s’®*’^)  be  obtained  at  90  degree  zenith 
viewing.  These  measurements  will  require  accurate  site 
measurements  of  atmospheric  data.  The  purpose  of  this 
experiment  would;  be  to  separate  errors  due  to  the  spherical 
refractive  geometry  and  the  transmission  modeling  at  the  high 
zenith  angles,  indicating  which  part  of  the  computer  model 
needs  revision. 

This  experiment  should  be  perforaed  under  controlled 
conditions  where  the  atmosphere  and  path  are  well  defined. 
Such  a  path  could'  be  achieved  =by  setting  up  a  path  between  two 
locations  and-  fixing  mirrors  at  each  location.  By  arranging 
the  viewing  path  such  that  it  was  bounced  back  and  forth 
between  the  buildings  a  series  of  long  paths  could  be  obtained 
(by  controlling  the  number  of  reflected  paths) .  Atmospheric 
data  could  be  collected  for  the  path  between  the  buildings, 
removing  the  other  source  of  non-model  error. 
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The  derived  atmospheric  transmissions  are  presented  in  this  chapter. 
The  original  solar  spectral  radiance  data,  collected  by  AFGL,  is  compiled 
in  Appendix  A. - 

The  following  table  summarizes  the  major  absorption  regions  and  the 
absorbing;  molecule: 


Wavelength 
Region  L) 

2.0 

■  2.75 

3.3 

4.3 
4.8 
5.5 


Primary 

Atmospheric 

Absorber 


H2O 

H2Q 

CH4 

CO2 

CO2.Q3 

H2Q 


Figures  4  through  8  are  the  atmosphenic  transmission  at  the  remaining 
zenith  angles.  The  same  general  features  described  above  are  evident  in 
these?  The  overall  transmissi bn  decreases  as  the  zenith  angle  increases. 
At  the  largest  zenith  angle  there  is  no  transmission  above  0.8  and  most 
bands  usually  considered  atmospheric  windows  show-  average  transmissions  of 
0.5  or  less. 

2.2  lComparison_bf  U0WTRAN6: to  Data 

The  models  were  compared  to  the  data  by  calculating  the  percent 
difference  between  the  data  and  nradel?  The  percent  difference  is  defined 
as: 


A  jj(X)  = 


T.(  X)  - 

c  m 

T  { A)  ' 

m 


where  t^(A)  is  the  calculated  transmission,  and  t^(A)  is  the  measured 
transmission.  A  positive  percent  difference  represents  an  overestimate  of 


the  atmospheric  transmission  by  the  model  and  a  negative  value  an 
underestimate^  If  %(^)  is  less  than;  5it,  t%{\)  is  undefined  as  it  is  felt 
that  there  is  insufficient  intensity  to  make  an  accurate  measurement  of  the 
transmission. 

L0WTRAN6;  accurately  models  the  atmospheric  transmission  at  the  smaller 
zenith  angles.  Figure  9  is  the  percent  difference  between  LOWTRANG  and  the 
data.  At  61.9®  and  65.5*  zenith,  the  percent  difference  is  close  to  zero 
at  most  wavelengths.  In  the  H2O  and  CO2  absorption  bands  there  is  a 
significant  deviation,  but  the  transmission  is  very  low  in  these  regions, 
causing  a  small  (  5%)  difference  in  transmission  to  be  a  large  percent 
difference.  As  the  zenith  angle  increases,  so  does  the  uncertainty. 

there  are  a  few  noticeable  features  in  the  percent  fdifference  at  the 
larger  zenith  angles.  A  large  uncertainty,  with  a  distinct  Structure, 
exists  in  the  2. 1-2.4  micron  region.  "This  absorption  is  due  to  CH4.  This 
absorption  band  is  entirely  missing  from  the  L0WtRAN6  data  base  (note  thH 
it  is  not  ap^parent  at  the  smaller  zenith  angles).  We  made  a  modification 
to  the  LOWTRANG  data  base  to  include  this  band;  it  is  discussed  in  Section 
5.  On  the  edge  of  the  2.8  micron  H2P  and  4.3  micron  CO2  absorption  bands, 
there  is  a  large  difference.  This  is  primarily  due  to  low  transmission  as 
mentioned  previously.  A  large  difference  also  exists  at  4.75  microns.  Two 
molecules  absorb  in  this  region;  O3  -  4.75-^4.8  microns  and  N2O  -  4i8-4.9 
microns.  If  the  O3  concentration  is  adjusted  to  obtain- better  agreement, 
the  agreement  in:  the  3.6  micron  03  absorptioh  band  is  lost,  N2O  is  part  of 
the  uniformly  mixed  gases,  so  its  concentration  cannot  be  adjusted 
individually.  This  example  indicates  a  basic  modeling  problem  in  LOWtRANG. 

2.3  Comparison  of  LOWTRANSB jModel  to  Data 

LOWTRAN5B  is  a  mid  resolution  atmospheric  transmission  model.  It  uses 
a  band  model  parameter  formulation  to  calculate  the  atmospheric 
transmission,  this  approach  more  accurately  takes  into  account  the  changes 
in  the  transmission  due  to  different  line  shapes  at  different  t^erUufes 
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and  pressures.  Consequently,  it  is  a  more  accurate  model.  It  also  has  a 
data  with  a  resolution  of  5  cm"l=  as  compared  to  the  20  cm~l  resolution  of 
tOWTRANB.  Due  to  this  increased  resolution,  the  calculations  of  L0WTRAN5B 
more  accurately  represent  the  structure  of  the  atmospheric  transmission, 
much  of  this  structure  is  not  present  in  the  L0WTRAN6  calculation.  The 
transmission  data:  display  much  finer  structure  than  the  L0WTRAN5B 
calculation  (on  the  order  of  0.01  cm"^)  and  FASCOOE  must  be  used  if 
modeling  of  this  fine  structure  is  necessary.  Note  that  there  is  more 
detail  in  the  percent  difference  for  L0WTRAN5B  and  FASCODE  than  for 
L0WTRAN6;  this  is  because  of  the  higher  resolution  of  the  claculatiohs. 

Figure  10  indicates  the  .percent  difference  between  L0WTRAN5B  and:  the 
data.  At  the  small  zenith  angles  L0WTRAN5B  does  a  fairly  good  job  of 
modeling  the  atmosphere;  essentially  equal  to  ‘LOWTRANSi  The  percent 
difference  at  the  61.7®  zenith  angle  is  approximately  0  at  most  of  the 
wavelengths,,  again,  excepting  the  areas  of  the -CO^  and  the  H2O  absorption 
bands  where  the  atmospheric  transmission  drops  quite  low.  It  is  quite  good 
at  the  65.5®  zenith  angle  with  small  spikes  being  observed  in  the  region  of 
the  CH4  .absorption  at  3.3  microns  and  the  beginnings  of  ah  inaccuracy  in 
modeling  the  ozone  at  4.75-4.8  microns.  This  apparent  modeling  deficiency 
is  of  the  same  nature  as  that  mentioned  for  tOWTRANB.  Attempting  to 
correct  the  ozone  content  to  make  it  more  accurately  model  the  4.75  micron 
region  decreases  the  ability  to  accurately  model  the  ozone  absorption  at 
3.6  microns.  At  a  zenith  angle  of  87.3®,  the  methane  absorption  at  3.3 
microns  is  more  evident  as  being  an  area  of  uncertainty.  There  is 
significant  structure  to  the  CH4  absorption  in  this  band.  L0WTRAN5B  does  a 
significantly  better  job  of  modeling  this  structure  than  tOWTRAN6  which  is 
inherent  in  the  designed  resolution  of  the  respective  models.  At  2.3 
microns,  the  sharp  positive  spike  is  due  to  methane  absorption  which  is 
completely  unmodeled  in  the  1.0WTRAN58  data  base.  This  was  mentioned  in 
L0WTRAN6  as  being  a  problem  and  the  problem  increases  with  increasing 
zenith.  At  a  90.3®  zenith  and  at  4.8  microns  the  inaccuracy  of  modeling 
the  CO2  absorption,  band  is  beginning  to  become  evident  as  a  positive 
percent  difference. 


ZenitH:  ■  91i.7° 


Zenith;  =  90>3° 


Zenith:  =  88.8° 


Zenith  =  87i.3° 


Zenith  =  65.5° 


Zenith  *  61.9° 
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The  investigation  was  undertaken  to  determine  the  potential  effects  of 
uncertainties  in  actual  C02  concentration  (relatively  minor),  the 
difference  between  the  band -model  codes  and  the  line-by-Tine  codes 
(significant)',  and  the  effects  of  convolving;  the  different  code  spectraT 
transmissivity  predictions  with  the  spectral  radiant  intensity  of  the 
plumes  The  investigation  has  been;  performed  for  a  number  of  scenarios  and 
the  complete  results  will  be  presented* in  a  report  which  is  currently  being 
prepared.  This  effort  has  been  exclusively  devoted  to  the  understanding  of 
the  atmospheric  transmission  effects  in  the  4.3-4. 6  micron  region. 

Near  the  end  of  the  computational  iphase  of  the  effort  we  realized  that 
the  optical  paths  for  each  of  the  target  altitudes  were  very  similar  to 
those  encountered  by  AFGL  during  spectral  radiant  intensity  measurements  of 
the  setting  sun  from  an  altitude  of  29,000-39,000  ft.  Therefore,  we 
compared  the  code  predictions  with  the  AFCL ‘measurements.  The  nomenclature 
and  conditions  are  indicated  in  Table  1. 


Table  1 

Aircraft 

Zenith 

AFGL  Run 

Altitude 

Anqle 

1 

39,000  ft 

61. 9 

5 

29,000  ft 

65.5 

27 

29,000  ft 

87.3 

Figure  2  indicates  the  comparison  between  the  actual  data  and  the 
predictions  employing  FASCOD2,  LOWTRANBB^  and  L0WTRAN6  for  Run  1  (similar 
to  a  target  at  30  kin  altitude  in  our  computational  scenario).  Clearly,  the 
L0WTRAN6  computations  do  not  describe  the  spectral  characteristics  of  the 
atmospheric  transmission  as  well  as  either  FASCOD2  Or  tOWTRANSB.  In  fact, 
L0WTRAN5B  gives  slightly  better  agreement  than  FASCOOZ. 

Figures  3  and  4  indicate  the  results  of  the  comparison  between  the 
predictions  and  data  for  Runs  5  and  27.  For  both  cases  L0WTRAN6  does  not 
provide  an  adequate  spectraT  description.  In  Run  5,  both  FASCOD2  and 
tOWTRANBB  provide  very  similar  results,  both  quite  good,  ^or  Run  27, 
FASC0D2  provides  much  better  agreement  with  data. 


A  further  test  was  conducted  by  convolving  the  actual  transmission  and 
the  computed  transmissions  with?  a  plume  spectrum  and  a  samp:le  sensor 
spectral  fiflter,  The  results  are  provided:  in  Figures  5^  6,  and  7  for  Runs 
1,  5,  and  27,  respectively^  These  results  echo  those  stated  in  the 
previous  paragraphs  -  namely^,  LOWTRAN6  is  not  to  :be  used,  L0V!TRAN5B  does 
better  ac  higher  target  attitudes,  FASC0D2  does  better  at  low  target 
altitudes. 


•  L0VVTRAN6 
llbWTRANSB'  i 

•  FASC0D2 

•  ACTUAL  DATA 


Figure  2.  Comparison'  of  'Atmospheric  'TransRii^ssion  Predictions  with  AFGL 
Data  Run  1  Zenith  Angle  of  6i.9  Degrees. 
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Figure  3.  Comparison  of  Atmospheric  Transmission  Predictions  with  AF6L 
'Data  Run  5  Zenith?  ^Angl  e>  of  '65.5  Degrees . 


•  L0WIf=l/!kN6 
■  L0WTRAN5B 
■FASC0D2 


fi;cUi!e>  4.  Comparison  of  Atmospheric  Transmission  Predictions  with  AFGL 
Data  Run  27  Zenith  Angle  of  87.3  Degrees. 


AFGL/prr  Data  Tape  Pornat  ^7  0 

Description  of  Kagnetlc  Tape  Fornat 

,  ibOD 

Tape  Is  lH  inch,  AOO  bpl,  9  track  and  consists  of  flles  ^separated  by  EOF 
narks,  with  two  EOF  ncrks  at  the  end  of  the  data. 


Each  file  consists  of  nultiple  records,  the  nuad>er  depending  on  the  type 
and  anount  of  data  In  .the  file.  Data  will  be  of  two  types;  spectral  and 
spatial.  Ihelr  respective  formats  will  be  discussed  separately* 

Spectral 

The  first  physical  record  In  a  spectral  data  file  will  be  a  180  byce 
ASCII  header  record.  The  ASCII  header  Is  broken  down  as  follows 't 


1  to  1'  bytes,  ir  bytes  r  A  nuitlpllcative  scale  factor  for  the  binary 

integer  data  array  in  E11.5  f/>nBat 


12  to  19  bytes,  8  bytes  -  The  time  or  frequency  spacing,  of  the  data 

5,rray  in  F8.5  fornat 

20  to  27  bytes,  8  bytes  **  For  frequency  data,  the  starting  frequency 

of  the  data  array  in  F8.2  fornat 


28  to  32  bytes, 


5  bytes  The  mmber  of  data  values  in  the  dsta  array 
in  15  forust 


33  to 

Ji. 


47  bytes. 


15  bytes  •  The  AFGL  nisslon  designator  in  the  ASCII 
fornat;  “MISSION*',  3A2 


48  to  57  bytes. 


1C  bytes  r  The  AFGL  archival  tape  designator  in  the 
/“Cn  format  ;  ~14re",^Ar / 

11  bytes  Tne  AFGL  archival  tape  rccoyd  nuab&r  or 

Qist^r  ftrvice  Nuid>er  (CSM)  in  the  ASCII 
format;  "RECORD  .  or-  “  CSN  -  '*  13,  If 
this  file  is  constructed  froa  more  than  one 
archival  record  the  record  wnber*vll]  be 
zero  and  the  source  records  description 
will  be  placed  in  the  last  108  bytes. 

4  b^cs  ~  ASCII  blanks 


73  to  bytes,  108  bytes  r  The  AFGL  archival;  record  descriptor  in 
it 3  ASCII  “A"  fornat  filled  with  blanks 

IN-  Y/O 

The  renal ning  records  in  the  file  are  binary  records  2880  bytes  long. 

They  consist  of  1440  16  bit  words  in  twos-conpllnent  fornat.  Each  word 
represents  one  data  point  in  the  dsta  array.  If  the  record  contains  less 
than' 1440  data  pdints^  It  is  zero  filled.  The  nu^er  of  data  records  in 
the  file  can  be  deternlned  by  dividing  the  nuiber  of  dsta  values  as 
determined  fron  bytes  28-32  in  the  header  record,  by  1440  and  founding 


to  the  next  larger  integer*  Ihe  value  of  each  data  point  can  be  determined 
by  multiplying  the  acaie  factor  from  header  bytea  1~’11»  by  the  aigned 
integer  value  in  the  data  record.  For  ;frequency  data,  the  frequency  of 
the  first  data  point  and  the  frequency  increment  per  data  point  aVe' 
contained  in  bytes  i2>27  of  the  header.  For  time  domain  data,  the 
interferometer  path  increment  is  contained  in  bytes  12-19  of  the  header* 

Norelco  iWlfc  foatial 


Ihe  first  physical  record  is  also  a  180  byte  ASCII  header  record.  Ihe 
header  is  broken  down  as  follows : 


Bytes 


Format 


1-11  \  El 1.5 


Description 


A  multlpli cativc  scale  factor  for  the  data  array 


12^16  F5^ 


Ihe  vertical  angular  spacing  of  the  data  points  within 
a  scan  in  milllradians/ 
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#  Iscans/reccrd 
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Appendix  3 
Run  01  Numeric  File 
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1  -  ■-■I .  i? 

?  cl  - -  5114911  7 

5-.  2?.?1 

3 ,  &7 1  cs  ?  5  7-  '  71  67  £  -  3  X'  v  - 

ivoi.AO 

X  .461  tli-ilio:  £36£- jX  07 

•  4.  -»  C  . 

5 .  j  6  3  1  04  266  ‘  6  6  E"  C  X”-  4 

1 905.33 

3 . 2  5  9  30  £  5  7  79  9  50  CE -0  0  07 

5r.24£4 

3. 0021(5  2S462  671 3  £-03(04 

1907.25 

1.231  41-9  9  4  3  9  9-3  6  8  £  -  0  0  0 7 

5.2431 

2. 9 76372 63 4734 2 9 E- 000 4 

1909.18 

■  3.303264S4799540E-0007 

5.2378 

3.0236^1 72886094E-0004 

1911.11 

8  .1 860201 94^99e76E-^0C07 

5,.  2326 

2 .9  9  2  23  S  4  4  3  3  95  5  8  E-  0  004 

1913.04 

3  .1 5:2  9  4  9  37  0(9  3  546  £  -  0  Q  O  7 

5..  2  27X5 

2  .98  015 0  0  5  09  61  £ 6  E-0 07  4 

\j\  y» 


.49v3rG=ji' 

.70=3  6:91511 
,?7i32?'?5:9; 
.336770  61  7 
.£,94  31145  2; 


rvlPc 

»  7  -1  ..  I 

5^763 
9  3  £31 
9^  67  S 
6  *  7  -  6 


«>  w  V 

0  DCs 
0CO6 
oc  3;d/ 

3  3  36 
o.*  ,■.> 


1  90  "12-''^^' 

i;2  91i9034;3r 
44  749  615  9  9 
3.63;947066r9952S5 
3..  74^38-93 5;20240E 
3.694433700X0153= 
3.31;9643  559  99  215  = 
3. 72  6  204  071199  792  = 
3.20051  S094179504E 


OCOo 
DCO' 
0  CO  6 
00C6 
Scot 

-OCX  a 
-OC06 


-OCC  6 
-0026 


7 : 0  c 
,9  3;?  2 
3  = 

■J  ?  T  ,* 
*■ 

a  t;/  ? 

a  315  A 
,  9  3:4  9 
,  =  342 


.  V  -  i  V,' 

.05:13 

.  =  3  0  5 
..  9  2  9  3 
.  9  29 1 


r.927C 
1;.  =  2:62 
1.9255 
1\.3243 
i.924l 
1.9234 


9.3  294=54227C;5  =  3^ 
9.8  4  4502  5655  9710; 

9.4  5X975  21  354  5  02 
r.X  2276011  26  9122 
1.G  5 55  3:  93  32X3  64 
?.3  74'4?60-5110c1 
5.5154  75X564^4  513 
9.791:9h536  665967 
1.065X62352=3196^ 
1V.0  2312M;6  574  9c5 
9.4X4^14*47190437 

9.406^51  3455:5102 

9.7  =  64649133=01  = 
9  . XI  6  ■’3=32  £33  5  9= 
•?.3  4X5Z3:3C4C;i,?03 
C.4  7311 '  =5i.5  32  92 

9.5  5  ■  V2  =  3027s72:? 
9. 301:39 '0  =  372  2  72 

9. 7  5  5:3033724:530^= 
1  .X 1 '3345  54746  54 
7.703:134X0  26:2728 
1  .03c2974742'^2566' 
5.4:3  654417  6W3  67 

3.6  =  999490933092 


■j  .■  -  4 
:*  I  *  ? 

.I  '  7 

>  ^  ^ 
A  T 

'■•r-  '■  7 


'X023 

.  n  n  t- 


=-0X23 
=-::c23 
c-CGG  2 
E-0-:O3 
£-0X02 
£-0X33 


Appendix  6 

Run  01  Interpolated  Data  File  (cm'^  units) 


1  «>  W  w  •-«/  W-  •/  ..  «»-•  .  »^.>r  - 

=  .  4211365  £-37 

■1  ;. r  f -» -  n  r*o  ‘^0  ''  “ 

1.C6:G‘'j1:  £*:^ 

i;i;c.:303003;!?cr^T 

3.1  31731:9  f-G 7 

1  ^1:5 

7.  466776^  £-27 

4  -  r\  -•  /»*--S  <■-  -| 

Im*.  ••»«  *W 

L, l-^llrll  --;7 

•i  “  "J  _  '  A  *■;  "^  p  ^  P  -T  ^ 

?  -  1  /  J  ?  A  ^  •• 

^  :y  *  'vn'  -^r  Art-  -% 

1.1 2:34323  :-:7 

1  "  J5  •  D  ’*'  ^'r'D 

7.5 '^93742  £-27 

4  i-'  r-  .A“'**-'^,*  rt/S-‘^rtr  ‘"A 

1  *.  *4  L'  •  .  w  O  Ur  ft.  •-•-  .f 

?.:5^65’2£-"7 

4  „  /=  r  A  '•-  “  rt’^  <-  O’.*'  rt 

1  V  *  -•  * 

Oafir^JI  £^:-7 

^  "  £  c  i  T  —r  - 

•  #ft#  a.  .  * 

4  V.  c  r,  -*  rt  r  r^  r  -  ,^  -  -  ^  -  - 

a.  •_.  *•=  1  .  • 

4  r  i  A  ;?  rt  rt  p  rt  -j  r  -  A 

1  M  ^  a.w  w  ««-•/-•  ^  w=r«  •  •--*  * 

t  a  5 c- 93  3 ?,3-  2^37 

4  ^-t  »,  -i  —  rk  ■•/p.’p  A  rt  rt 

t  M>  •-«.*.  ft^  «*'  »*-"•>'  •» 

7  i -23  62:522  £-7;  7 

>  <n  •*-*. 

1  C.r  *  m-\J  ••  ‘J-J  \J  ^  •J  mf 

6.  5:1:561  £^27 

4  r  T  S’  A  -  A.  f-  n  A  «  Tt  p  r  p 

7  ^  '  5  '  7  r  -  '.7 

■*.••■,  'I  ''.''!"  "  "if '''.f:' T 

*  ,J  •  W  w  V  w  W  ^  w  W  V. -W  V 

is  5^^1"}2^'7 

■1  s  i  nriiif;  f' '■.''r. 

1.C*'37  5  55  1-'^ 

1  .'  *v,a  -J  J  •- J  o..- v  ;  -; . 

1.9f^1  6"57-~7 

1  S5  5  a  C'OCD  C003  DCGO 

£.12627  65  £-■::’ 

1  PDOaCGOCGODO^CGC 

^.41  971  £51^27 

1  905  a  1003  00002,0  :■ 

c  r3C*  =  7i,  ri"7 

IPTOaCGOlGOOlGlOO 

2. 2657301 £-27 

1  91  5  a  00000003 C GOO 

?.34;£74i3r^G7 

152.C.OOGGCDD3CG:00 

7i '591  l  ea  £-27 

T925.OO0000a0OG0O 

8^1900464 £-07 

1930.000000000000 

8.7557336 E-07 

1935.0000.00000000 

7i 9999239 Et07 

r940.o;ooo.oooooc::00 

S>  9.674813  £^07 

1945.O00OG0G0GGOC. 

8 i 4756374  £?G? 

1?5C.O0CO0020GO;CO 

8.5933291  S^27 

1955.000000000000 

Si  8629764 £r37 

196O.G3C3000CO00Q 

gv95e92  54£^07 

1965.000000000030 

9.0e'496^cE^07 

1910.030333000003 

£.9OC1099?40? 

1975.000003300003 

9^376331 6 E-07 

1930.000000000003 

9.1 945527 £t37 

1935.600000000000 

9i37973T4E^07 

1990.000000000600. 

9. 13 0006? E-07 

1995.000000600000 

8. 9944503 E^O 7 

2000.630300000003 

9. 2475393  £^07 

2005. 600300000003 

9.4966612 E-07 

201G.6003006oc600 

9V43:9l953E-07 

2015.000000000000 

9. 6141 787  E-07 

2020.000600600000 

9.33.10496=^07 

2025.000000000600 

9i 471 2962 E-07 

2030.000000000000 

9.5223821 E^O? 

2035.020000600000 

9.  6291  501  £^07 

2040.60G30330GGC3 

9..  41  3  53  55  =-37 

2045.30C00360C0G0 

9^  36^93755  E-G7 

2C5C.60GOC36o:-3G3 

4.22G7915E-:-7 

Z 355.303605312232 

S.  7430430  =^''7 

a^  f*  ^  -»  17  p  A  A  p,  A  A  C,  A 

^  AACp.  T 

^  -•  -  .  .#  »  ■#  -a  .  - 

1 0i  3  •  ^  i  loi  13  T  3  r 

1.  5£c51  f2  2-''7 

.•  W  ~  W.*  ft#  W  W^^ai#wWW‘» 

2u75.330600600GGD 

5.  5976313 --DT 

2080. 000600600000 

6. 3176407 =-D7 

2055.000600600600 

8. 4590529 E-07 

2090.600600000600 

6.3752935  £-37 

2095.000000000000 

8.1523571  E-07 

il>05. 
3r-io. 
5^  1  5 . 
f.1  2  j  . 


coo30')C3;oc::i 

0  303C'00::03D-D 


;•  1 


31  3  5  i 


;  1  i  C . 
3T45. 
rt;;. 
5^  3: . 
Z‘  -1  0  'J  • 

5:1  65. 
31 T  u . 
5173. 
513  C . 
3135^ 


31^0:. 


31?  3:. 


00^:Ci00D:G0o3 
C100000?1}OOG 
C7G30003;03C§ 
3:OGOOCij232} 
3J3  3^^30  30  03) 
3  3G3(y00D:33jG 
D3^303CQ^bC33 
3  3:3  3  03)  3  0:3  0G3 

53330:303:0^00 

0  73303030C&3. 
030333330033 
j0;0  3  0300:000  3 

w'  w  ^J  w  «/.U  ^  «« 

GIS 3GD00iCC5 

'  1  n  T'  ^  A  r.  -» 

«  *<  w  ^  *^-V*  ^  v  ' 

333300000033 
r  ^  ^  n  nr  .n  n  n  n  n 

••  w*  ^  ^  w-w  w  U  «.■ 

300300030000 


7.31  351  ;2‘=- 
3.^413^0^6  0“ 
3.G533t.:S  “- 
3.23  6:3<>40r“ 
3.05?4305  0- 
3.:5C04;ti.:- 
3.53  01 '73  > 


-Ji 

:t 


“3-! 


i .  i  y  i  c  y  *1  r  “ 

2.40  2121:^  =, 

7.31  ^OO-iO' 


•  --«  /• 
.  r  f 


.73131  31  2-' 


3.36:7556“ 
3.  9  6  5,3  56?  Z 
3.::  0  73  71  5“ 
3. 54337^1  E 
3.57235573 
3..c?63462r 
3.3 014  40c,  r 


-  n  y 


Appendix  7 

Run  01  Interpolated  Data  FMe  (micron  units) 


5^5556 
5.5402 
5.5249 
5.5096 
5.4945 
5.4795 
5.4645 
5.4496 
5.4346 
5.4201 
5.4054 
5.-3908 
5.3763 
5.3619 
5.3476 
5.3333 
5.3191 
5.3050 
5.2910 
5.2770 
5.2632 
5.2493i 
5^2356= 
5.3219; 
5.2083^ 
5.1948 
5.3813^ 
5.1680 
5.1546 
5.1414 
5il282 
5.1151 
•5^1026 
5.1891 
5.0761 
5.1633 
5.1505 
5.1378= 
5.0251 
5.0125* 
5.1000^ 
4^1875= 
4.1751 
4.16i28i  i 
4.9505 
4.1383 
4*9261 
4^1140> 
4.1021' 
4.19O0 
4.58780 
4.8662 
4*8544 
4.8426 
4.130  9; 
4^1193 
4*8071 
4.7962 
4.1847 
4.1733 


2 . 736 14 7 4  61 6  8 16  0  E -1 0  0  4 
3 .44  623 1844388  41 E- 0  0  0  4 
2i64916325576792E-0004 
2.4  5294805028884  E-10  04 
2.9  2718332189246  E -0  C  0  4 
2i37297160S58203E-0004 
2.37804  637632566E-6i 04 

2  *5  3  8  3  7 4  911 3  0  94  8  E - 00  0  4 
2*  48  2  21 2 1 8797973E-0  004 
2.06186338552961E-0004 
2il3311328229577E-1004 
2.17  6  4:1 420239945E  -11 0  4 
2.26670505224735E-0004 
2.44081 5239 13753  E-1 0  04 
2^ 2951 03978 04297E^6 004 
2,i  7  5  424  7  6  21 7 1 34  6  E-  6  0  0  4 
2.773438343695636-1004 
2*861036075882956-1004 
2.83379016730922E-0004 
2. 9748  41 50249496E-1 004 
3.03151975418636E-6o04 
3.00043975  8  347 18  E-Ol 04 
3^0  0861 9045852266-01 04 
3*053677361579696-0004 
2*900i^70464422'99Erl004 
3*0276411f900187e-0004 
3.264119676091866-^6014 
2 i9  8761158019 19  0  E- 6l 04 
3.36630280115821E-6004 
3*19811225904854E-6o04 
3 . 3  73  6 17  3954449  3  6-0004 
3*3787  8  8lt687987E-6  004 
3.43211512344960E-^6i04 
3*498984423901956-0004 
3*445  277  04222219  6-0004 
3^48696215509436-1104 
3*595647149969856-0004 
3*6865 1 5  8  3329  8  8  2  6-1 0  04 
3.606489673377846-6004 
3*570853653478426-6004 
3.687773381249846-6004 

3  *  8  08161138407766-1 004 
3 - 8  0404  29  63778  896-1004 
3*  893886584174306-0104 
3*798017121012806-6104 
3. 874234528860046-6004 
3*914824294848356-0014 
3*977850054264216-6004 
3*907936371083666-6004 
3.908717323559696-6104 
3^865586316251^86-0104 
3*683183312785946*6004 
3*685115066514536-0004 
3 i 6  9  3 3  3  9  478  9  6  067  6-01 0  4 
3  *  648  9760 1826  45  36-10  04 
3*  8 187 4840  78  9 99  7  6 - 01 0  4 
2.726693726260446*6004 
3.668535072951516-1004 
3*649661950486(666-0004 
3.569563454237656-1004 


4.7619 
4.7i06 
4.T393 
4.7281 
4.7170 
4.7059 
4.6943 
4.6838 
4.6729 
4.6620 
4.6M2 
4.6^04 
4.6296 
4.6189 
4.6083 
4.5977 
4.5672 
4.5:767 
4.5662 
4.5558 
4.5455 
4.5351 
4.5i49 
4.5147 
4.5045 
4.4944 
4.4843 
4.4743 
4.4:643 
4.4'543 
4.4444 
4*^346 
4.4’248 
4.4150 
4. 40 53 
4.3956 
4.3860 
4.3764 
4.3668 
4.3573 
4.3478 
4.3384 
4.3290 
4.il97 
4.3103 
4.3011 
4.2918 
4.2827 
4.2735 
4.2644 
4.2353 
4.2463 
4.2373 
4.2283 
4.2194 
4.2105 
4.2017 
4.1929 
4.1841 
4.1754 


3V63832675388309E-0004 
3 .71 5  9  08  327257  54  E-0  0 0  4 
3.56686i42099326E-0004 
4.14258633242959E-o604 
3 . 50  73  00  997  274 14  E-0  00  4 
3.40159165542004E-0004 
3. 73  4  8042  2  829  6  97  E- 0  00  4 
4.46410314  08  28  73  E- 0  004 
4 .17 6 976  8  4  24  9  9 8  9  E-  0  00  4 
4 .91 5 4  8  8  3  275  4  5  9  6  E-0  6  0  4 
4  .7.7  8  6  49  2  328  46 10  E-  0  0  04 
4 . 7  6  5 1 3  0  8  2  8  2  28 16  E-  0  6  04 
4 .76 48487  817  9  3  06  E - 0  0  0  4 
4.73659883li6095E-0O04 
4.74536077258225E-0004 
4.7S119391791345E-0004 
4.608532929:63304E“0604 
4.51580107350580E-0004 
4 .19  85047  281 567  5  E-0  0  0  4 
4 .0  3 6  3  42  0 18  2  5  2  09  E - 00  0  4 
3 .81 772197367347  E-0  004 
3.S0159053903010E-0004 
3.29570843689488E-0O04 
3 . 03273  51 1154778  e^OO  04 
3.07358216494968E-0004 
3.95 8 15 17  0  3 6  824  lE-0  6  0 4 
3. 772968867871218-0  004 
2.504996937i249iE-o604 
2 . 1234  3 106723  800  E-0  004 
2.i8073896881554E-0004 
2.36823638584793E-0004 
1.81194135368212E-6o04 
1.38409662550742E-0004 
9.32264203673538E-0005 
4.33904186120659E-0005 
1 . 4 7167 6  6 1 1 7  6  0  5 4E - 0  0  0  5 
1.53723804625063E-0005 
1.0  0  357 1744  9  6191E-0  0  04 
2.67344146779869E-0005 
1.0721968951313t3E-0005 
4.70322125262213E-0006 
4.35103637085738E-0006 
l.:6556U89660400E-0005 
.7.52567O3O28OO66E-0OO6 
9.56959207591701E-0006 
5.08753563016324E-0008 
3.51333858080935E-6006 
44885836202310i3E-0007 
1 .223134339348528-0006 
4.94436474970938E-0006 
1.93371 3744066428-000  6 
4  ^87  6  07  67  5  U  9 113  E- 00  0  6 
4.21883355385816E-0006 
1  i^^8  3  646  5  6794  694  9  E-  6  0  0  5 
4. 4182517  3344568  E-0  0  0  5 
4 .2933842  8 04 1270  E-0  005 
3.589922082292i4E-6005 
1^56728949333829E-b004 
545528i668577656E-0004 
7il2293226338367E-0004 


4.1667 

4.1580 

4.1494 

4.1406 

4.1322 

4.1237 

4.1t52 

4.10^8 

4.0984 

4.0900 

4.0816 

4.0733 

4.0650 

4.0566 

4.0486 

4.0404 

4.0323 

4.0241 

4.0161 

4.0080 

4.0000 

3.9920 

3.9841 

3.9761 

3.9683 

■3.9604 

3.9516 

3.9448 

3. 93 JO 

3.9293 

3.9216 

3.9139 

3.9062 

3.8986 

3.89il 

3.8835 ‘ 

3.8760 

3.8685 

3.8610 

3.8536 

3.8462 

3.8388 

3.8314 

3.8241 

3.8148 

3.8095 

3.8023 

3.73h 

3.78/79 

3.7807 

3.7736 

3.7665 

3.7594 

3.7323 

3.7453 

3.7383 

3.7313 

3.7244 

3.7t75 

3.7106 


7 . 53  8  9  606 8758309 E-0  0  04 
7. 60675199529808E-0004 
7. 50687935606287E-0004 

7 . 3  77 9  6  61 0  5  3  3  0  8 1 E  - 0  0  04 

7 . 4  2  5  94  5  2  5  26  08  2  0  E-0  0  04 
7.49  47658  0351072  E-0  004 
7.  67927269716751E-0004 
7.73746061 244474E-0  004 
7.76324118550775E-0004 
7. 986861511 96826E-0  004 
7 . 943  29194 15  81 17  E- 0  0  0  4 
8 . 1 04  22241 8  01 9  4  5  E- 0  0  04 
8 . 121 2  8707  8  335  37E-0004 
8.16402564537633E-0064 
7. 92166262735705E-0064 
8 . 312  8597045956 IE-0004 
8.502384607513856-0004 
8 . 660022645  05  2  83  E-0  004 
8.67871992105584E-00d4 
8.699210108131796-0004 
8.84705855590973E-0004 
8.83496966828545E-00d4 
8 .863  2 0273  03958  7 E- 0  00 4 
9.029l4931622867E-00d4 
9 . 16  9  2  9  6447  0  8  0 2  0  E-0  0  04 
9.250484215108661-0604 
9.24525938814647£-00d4 
9.3^9271683888866-0064 
9.27611049839072E-0004 
9.18760874436764E-Oo64 
9.230074194892526-0004 
9.2472265859:9997E-0004 
9 .0420  0  8459  303471-0  00  4 
9 . 59 154  507  83i23  8E-0  0  04 
9. 333 8 2 268 978 10 06-0004 
9.38486686I35498E-O064 
9.4655943638O667E-OO04 
9.579OO6879532O6E-O064 

9. 84241481440051E-0004 
1.030232923049826-0003 
1.0402li23595466E-0003 
1.0491 18  397  2920  5  E-0  06  3 
1. 060478627678666-0003 
1. 064467047678266-0003 
1 . 0771 71320 S59  23  6-0  00  3 
1.089660714826676-0003 
1.094353109321496-0003 
1.09870467944795E-0003 
1.103O21879528O4E-OO03 
1 . 11 50977166742  8E-0  003 
1.102090592153716-0003 
1.108426830358066-0003 
l.l4908684660620E-00d3 
1 . 1 4  0 177966  89  5  81E-0  00  3 
1.128285795688916-0003 
1.122910186714076-0003 
1 . 11 63  9  554  09 1449  E-0  00  3 
l.m605828l79446-0003 
1. 165154503166126-0003 
1. 15 1679243978  526- 0  003 


3.7b  3  7  1. 17  6 1047699  83  87  E-f  0  0  03 
3^6969  1. 169^  2133l95629E-000t3 

3.6905  i.  2  03  3  47  0 9  5 1 29  8 1  Et^O  0  0  3 
3.6832  1.2003056 102  24^1 1 E-O  003 
3  4  6765  1 . 185  2  26  35  56 141S9  E^OO  0  3 
3 .66  97  1 . 2  0  0  2  3  6  5  01 9  9  5  8  0  E^O  0  0  3 
3 .663  0  i . 20  60718 1 73  9  0  9  6  E^  00  03 
3.6563  .  i.23925615571707E-0003 
3.6496  1 . 256  5  63 13836409  E-^O  0  03. 
3 .6430  1 . 224449  062453  56E-0  003 
3.63  64  i . 2  52  287  0  4  8  0  9  64  3  E- 0  0  0  3 
34  6298  1.26019902873331E-.0003 

3 . 6232  i . 248  0  424 1 38  7544  E^O  0  03 
3.6166  1 . 2  8  3  7  7  9  2  9 1 2  3  8  5  7  E - 0  0  0  3 

3 . 61 0 i  i . 3  0  0  2  40  3  079  0743  E-0  0  0  3 

3.6036  i.25569003167492E^0003 

3.5971  1 . 3  2  9  8  6^  0  82 1 7  3  06  E-0  0  0  3 
3 .5  9  07  1 . 282 1 18  2  224  23  63  E-0  003 
3 .5  842  i . 3  0  8  2  3  5  29  04  074  7  EiO  0  0  3 

3 .5 77  6  1 . 2  8  3 6  54 4  7  7  8  73 6  0  E- 0  0  0  3 

3.57 14  1 . 293  2127  877  8  258  E-OO  03 
3.5651  1 . 31 1 3  647  080444  9  E-0  0  0  3 
3.5587  1.34953725056519E-OOO3 

3.5524  i  .  36  647 4  4  5  7  4  8  57  2  E-4 0  0  03 
3.5461  1.36505927i28283E^0003 

3.5398  i.35283223351124E-0003 

3.5336  i. 35441114375112E-0003 
3.5273  i.41346653065888E70003 

3.5211  1.40467630895458E?-000:3 

3.5149  1 . 4 4  5 4  8  57  5  5  3  2  61 0  E-0  0  0  3 
3  .15  0  88  1 . 41 65 17  2749  0  8 1 0  E-0  0  0  3 

3.2  026  1.4 1 8411278 166  33  Er>0  0  03 

3.4961  i.45294S56623377E40003 

3.4904  1.43729388464919E-0003 

3.4843  1.43691247621857Er0003 

3.. 4783  1.4458676109  008  2  E4 0  0 0  3 
3.4722  1.48006680787238 E^0003 
3.4662  1.4§576:230487305E40003 

3.4602  1.49012823719197E40003 

3.4542  1.48O6445O048224E?OOO3 

3.4483  i.48215634633253E’^0003 

3.4423  1 . 48  066  8  88  93  83  81 E40  003 
3.4364  1.45452161775594E40003 

3  .4 30  5  1  .i572 30515354 10 1E40  0  03 
3.4247  1.  3  3 13 16  3  377 4  3 97  EtO  00  3 
3 .41 8  8  1. 5807  646  073  9132  E^O  003 
3.41 3  6  1 . 345 482  8  8  5  3 5  0  0  8  E40  0  03 
3.4072  1.56433906605180E?0003 

3.4014  1. 34643992272387E-0003 

3.^  3956  1 . 56428068794412E='0  003 
3V3898  1.32967360295311 8^0003 

3 .3  8  41  1.598 199 17 1 7648  5  E-0  0  03 
3.3784  1.3208095 71 60 4:19  Er 0  0 0  3 

3 .3  727  1 . 64  23  75  8  2  5  5  23  9  6  E- 0  0  0  3 

3.3  670  1 . 43  5970  6  0  875 14  5  E-0  0  0  3 
3.3613  1.67921873659616E-0003 

3.3557  1.4753d014i60502E4O003 

3.350 1  1.  6731 49 266  9 93  69  E^O  003 
3.3445  1.48221153095562E-0003 

3.3389  .1.65882450743915E-0003 


3.3333 
3.  32^8 
3.32l3 
3.3167 
3.3^13 
3.3058 
3.30b3 
3.2949 
3.2895 
3.2841 
3.2787 
3.2733 
3.2680 
3.2626 
3.2573 
3.2520 
3.2486 
3.2415 
3.2362 
3.2310 
3.2258 
3.2206 

3.21:54 

3.21*63 

3.2051 

3.2000 

3.1949 

3.1858 

3.1847 

3.1757 

3.1746 

3.1656 

3.1646 

3.1596 

3.1546 

3.1496 

3.1447 

3.1357 

3.1348 

3.1299 

3.1250 

3.1201 

3.1153 

3.1104 

3.1056 

3.1008 

3.0960 

3-0512 

3.0864 

3.0817 

3.0769 

3.0722 

3.0675 

3.0628 

3.0581 

3.0534 

3.0488 

3.0441 

3.0355 

3-0349 


1.56501682394605E-0003 
1.761755293316016-0003 
1.550231933066476-0003 
1 . 1 6  8 1 764  5  82  97196-000  3 
5.35507250057243E-0004 
1.639361638128286-0603 
1.436103881860406-0003 
1 •72 84  27  8  51 475  05  6-0  0  0  3 
1 .73  4  357  5  8  623495  6-000  3 
1  •  8 4  5  0  2 8  2  5  6 3  8  8  6  6  6 - 0  6  0  3 
1 .59  8  8 119755  84986-0003 
1.687211488469846-0003 
1 .52  3365500901266-000  3 

1 .84  8653124710  61 6-00  0  3 
1.607215690331296-0003 
1.87431859  005645  6-00  0  3 
1.66026232^914726-0003 
1.82756272  00  4  5  25  E -00  0  3 
1 . 75  7  6 2 6  7  6  44 5  3  2 4  6 - 6  0  0  3 

1 . 79  4785  403  589576 -00  0 3 

1 .84  0 1 25  9  887  7  2  81 E-QTO  03 
1 .75  6  0  7327  67 1635  6-00  0  3 

1 .85  5814937478596-0003 
1.740402015252  946-00  0  3 
1.981290117.949906-0003 

1 . 80  36328  05025056-0003 
2.007186936545226-0003 
1.525794283655696-0003 
1. 569  5  83  5  5  5347  04  6-000  3 
1.95935'304157666E-6003 
1 • 951 628748402976-0003 
2.126496356378746-0003 
2 . 6747623378  5542  6-00  0  3 
2.09650687605336E-Q003 
2. 118358914078306-0003 
2.062228524078346-6003 
2.047746686866256-0003 
2.139586259573896-0003 
2.i343<491ll466ilE-0003 
2.154618526894766-0003 
2.108520057483086-0003 
2.160650384143456-0003 
2.1^5109673364616-6003 
2p 115613  3106  57746-0003 
2.156460860407366-0003 
2.194258038617396-0003 
2. 206887132619076-6003 
2.187070386622636-0003 
2.214179443220146-6003 
2.201188197524576-0003 
2.239264386868636-0003 
2.265030551036776-6003 
2.296272905802256-0003 
2.266171992321386-0003 
2.302766669263436-0003 
2.29177649982049E-6o03= 
2.369969649286536-0003 
2.336028571314056-0003: 
2.34  0952012537436-00  0  3 
2.^99665181855536-0003 


3.0303= 
3.0257 
3.6211 
3.0166 
3.0120 
3.0075 
3.  6 030. 
2.9985  . 
2.9940 
2.9895 
2.9851 
2.9806 
2.9762 
2.9718: 
2.9674 
2.9630: 
2.958  6. 
2.9542 
2.9499; 
2.9455 
2.9412 
2.9369 
2.9326 
2.9283 
2.9240 
2.9197 
2.9155 
2.9112 
2.9070 
2.9028 
2.8986 
2.8944 
2.8901 
2.8860 
2.8818 
2.8777 
2.8736 
2.8694 
2.8653 
2.8612 
2.8571 
2.8531 
2.8490 
2.8456 
2.8409 
2.8369 
2.8329 
2.8289 
2.:8249 
2.8209 
2.J169 
2.8129 
2 i 8090 
2.80  50 
2/801;l 
2.^797:2 
2.7933 
2.7894 
2^7855 
2i7816 


2. 3928  2  3062194796-0003 
Zi  47=83474700  95346-6003 
2 . 4  5  6  0  4 4^  8  30  5 15  3  E-6  0  03 
2 . 4  447  52  6 159  89646  -0  0  03 
2. 48  67  87473127  81 6-0003 
2.49445479162347  6-0  0  0  3 
2.5i'199S347315836-0003 
2.529854953419176-6003 
2.5196^3832446436^0003 
2.53912l9l55643eE-!0003 
2 . 5487493  76626  71 6-0003 
2.569325025652186-0003 
2 . 6027 1.0  5  0  76  6  292  6-0  003 
2.57^0422418262556-0003 
2.602476834276506-0003 
2 . 640535031556376 -0  0  03 
2 .7124 183 134  041 9  6-0  003 
2 . 6229 1 3  074  52311 6-0003 
2 . 672 1 9  7  2  9  4366196-0  003 
2 .678 174  235661866-0003 
2.733370048630956-0003 
2.7458 03685S06596?0003 
2.764395706428506-0003 
2 . 7433  3  8  4693  296 1640  0  0  3 
2.7452  74 35676930  640  0  0  3 
2.733695  5767403  0  640  0  63 
2.778151386063856-0003 
2.86590086316224640003 
2.853043664963206-0003 
2.735956585545776-0063 
2. 68163698644136646063 
2.788933996651976-0003 
2 .7721 942771 9  544640  003 
2.678332010077616-0003 
2.76987503380127640063 
2-76899945251912640003 
2 . 7749 61 6 01 95412  640  003 
2.65648592446510640003 
2.6098  937909 1 114  6?0  0  03 
2.61719581377231640003 
2.643379641437^35640003 
2.580341777633786^0003 
2.676689918043716-0003 
2.>58462045177387e40003 
2. 663912712812086^0063 
2.5370506  076  8539640003 
2. 66335561657090640003 
2.62407572072476640003 
2.50466133728366640003 
2 . 2844427  0  04  85566^0  003 
2.246417387734576-0063 
2.150962220746076^0003 
1.969325882864586-0063 
1.72  99  6  58  45225  5  9  640  003 
1.3  8376  2130156906-00  03 
1.5 344779  45  95955 1-0  003 
1.24667493604683640003 
9 . 418432968244526-0 0  54 
4.697760989729186-0064 
2.193147275528776-0064 


2.777 B 
2.7739 
2  •7701 
2.7663 
2.7624 
2.7586 
2.7548 
2.7510. 
2.7473 
2.7435 
2.7397 
2.7360 
2.7322 
2.7285 
2.T248 
2.7211 
2.7174 
2.7137 
2.7100 
2.7064 
2.7027 
2.6991 
2.6954 
2.6918 
2.6882 
2.6846 
2.6810 
2.6774 
2.6738 
2^6702 
2. 6667 
2.6631 
2.6’596 
2.6560 
2.6525 
2.6490 
2.6455 
2.6420 
2-6385 
2.6350 
2.6316 
2.4281 
2.6247 
2.6212 
2.6178 
2.6144 
2.  6110 
2.6076 
2.6042 
2.6006 
2.5974 
2.5940 
2.5907 
2.^873 
2.5840 
2.5806 
2.5773 
2.5740 
2.5707 
2.5674 


6.03134213842393E-0005 

4.13225775814374E-0005 

1.2924i61551S197E-0004 

7.59638739877389E-0004 

2 .5626  07945408  86  E-O  0  04 

1 . 06097976647179E-6o04 

1.58  780  6  623  2  266 6  E-0  0  04 

4.49680666896324E-0004 

1 .1 6 4 0  6 0  5  71 2  6 63  9 E - 0  0 0  3 

1 .8448085 1 3078  jB9E-6  003 

1.65767362931923E-0003 

2.4473465789  292  0  E-0  0  0  3 

2 . 61 5749  7  51 0  75  92  E-0  003 

2.74006808277960E-0003 

2 . 0523459  848  256  2E-0  0  0  3 

1.66919837964485E-0003 

1. 286867944 84700E-0 003 

1.19258003825884E-0003 

4.79136286458  01 5  E -0  0  0 4 

1 .64  828  2  6  57  3849  3  E-0  0  04 

2 .068  373402 1 83  69E-0  00  5> 

1.86588696530210E-0005 

l.iO  1 5  52  5  8  65  8  9948  E-0  0  04 

2.98056959638249E-0004 

2.64689400269489E-0004 

1.1468733 9.87 61 8  5  E-0  00  4 

1.033248 04205163E-0004 

1.028651192231.94E-0004 

3.54206086789155E40004 

1.11264943407541E-0003 

1.91767893142814E46003 

2 .5  8  79  7  5  0 5  0 1 4 03  0  E-0  0  0  3 

3.356 148  3  2292563  E-0  0  03 

3.75466734150019E-^0003 

3.46272  550475391 E-0  003 

3.83613833053T86E-0003 

3.45577295354005E-0003 

3.5963978145190*E-000^ 

3.68914351729899E^0003 

3 .940973 247  05456 E^O  003 

3^18862644547124E?^0003 

2.847133lil57248E-0003 

2 .9  0  5692297  24724  E^O  00  3 

3 .693369696868  84  E-0003 

2^92586322503041E^0003 

2.982i945e500881E-0003 

3. 5443266f297289E-0003 

3.40847788800502E^0001 

2.80969300582967E?0003 

3.34567843261269E40003 

3.8686228i054601E^0003 

2i 197631 13382854E^0003 

3.6401 3966240i5lE^0003 

3.31172606729i69E-0003 

3v396i 69  0141524 lE-O  003 

35503i2789908358E-0003 

4.08230696074696E?0003 

3i717395l8972447E-0003 

3s99253165979019E-000i 

3573433040957494E^0003 


2.56i41 

1.5608 

2.5575 

2. 5543 

2.5510 

2.5478 

2.5445 

2.54i3. 

2.5381 

2.5349 

2.5316 

2.5284 

2.5253 

2.5221 

2.5189 

2.5157 

2.5126 

2.5094 

2.5063 

2.5031 

2.5000 

2.4969 

2.4938 

2.4907 

2.4876 

2.4845 

2.4814 

2.4783 

2.4752 

2.4722 

2.4691 

2.4661 

2. 4651 

2.4600 

2.4570 

2.4540 

2.4510 

2.4480 

2.4450 

2.4420 

2.4390 

2.4361 

2.4331 

2.4301 

2.4272 

2.4242 

2.4213 

2.4184 

2.4155 

2.4125 

2.4096 

2.4067 

2.4038 

2.4010 

2.3981 

2.3952 

2.3923 

2.3895 

2*3866 

2.3838 


3.848  37  8 1 3 1 29  8  35  E-0  0  0  3 
3 . 279  8  8 57  3 191 0  Oil- 0  003 
4 . 15849  60 1 3  2 11 28  E-s-O  0  0  3 
4 . 28602076863882E~0003 
3. 930161237  88215eV0  003 
4.23077073516254E-O0b3 
4 . 156  5  250  4  5  58  039E-0  0  0  3 
3. 9744879 4467433E-0003 
4.48  495796372  5  e3E-0  0  0  3 
4.305040078122606-0003 
4.25279504229081E-0003 
4. 434250219397036-0003 
4.708910461353366-0 003 
4.565374392140306-0003 
4 . 6 1 3  5  2  37  26  0  21 646-0  6  0  3 
4.32835989066405E-0003 
4. 597398375  992166-0  003 
4.528015411722216-0(103 
4.779245310338356-0003 
4.647805199383016-0003 
4.645394948148866-0003 
4.624344888540586-0003 
4 . 64321083045  8  606-0  00  3 
4.643992953717206-0003 
4 . 87  6080670619846-0003 
4 .  667  8  73  28  7  6  0 8  8  7  6-0  00 3 
4.6211840734:99936-0003 
4.577880072801586-0003 
4 .9 3991 54 5  2  to 94  416-0  003 
4.7343li298{)7267E-0003 
5.015760857922666-0003 
4.99761739471438.6-0063 
5 . 0 129  9  f 1 424  5  8 1 6 6-0  60 3 
4.945172816690056-0003 
4.86516  6523472416-0  00  3 
5.0  3  024777651 0  076-0  00  3 
5.268931949461266-0003 
5.228885996238836-0003 
5.069909329662696-0003 
5.24996625121498E-0003 
5 . 1943  0  02466  99  876-0003 
5.15273701584817E-0603 
5.11707632864506e-0603 
5.178682855138566-0003 
5.2i085724179S26E-0603 
5. 062977559603616-0003 
5.361947251351756-0003 
5.335461750512146-0003 
5.218326556700966-0003 
5.420597797605666-0003 
5.4210l780325055e-0603 
5.337595964135746-0003 
5.44242245467785E-0603 
5.517697251207216-0603 
5.7  50  2  33  275  350076-0003 
5.3866866808  2920  E-OO  0  3 
5.825307032708346-0003 
5.364607337476466-0003 
5.57448012l4530lE-0003 
5.499561718679056-0003 


2.3810  5.33751377479064E-0003 
2 .37  81  5^63521977 76 26 5 6  E-O  00  3 
2.3753  5. 69368605069556E-0003 
2.3725  5.936162058873i5E-0003 
2.3697  6.95870062290038E-0003 
2. 3669  5.  73951966613606E>0003 
2.3661  5 . 86269  88 11 6  05  37E-00  03 
2.3613  .5.86763318571666E-0003 
2.3585  6. 08906881 375596E>0003 
2. 3557  5. 70519956965576E-0003 
2.3529  5.77011666358888E-0003 
2.3502  5.760811666110586-0003 
2.3676  5 . 7 1 5 65 6  2  2 6  3  72  25  E - 0  0  0  3 
2 .3667  5 . 66 16 3  9  2  33 1 226  8  E-O  0  0  3 
2.3619  5.96032369687962E-0003 
2.3392  6. 048009060 00236 E-0003 
2.3366  6.10597916593569E-0003 
2.3337  6. 02662975851587E-0003 
2.3310  6. 1352767161 0115E>0003 
2.3283  6.08069912302988E-0003 
2.3256  6.10657762025618E-0003 
2.3229  6. 19961138376713E-0003 
2.3202  6.13116323768931E-0003 
2.3175  5.96375066057318E-0003 
2* 3168  6.03bl053665606lE-0003 
2.3121  6. 196967J6236696E-b003 
2.3095  6.16123585675865E-0003 
2.306  8  6 . 61166663761665  E-0003 
2.30  61  5.88  821357 127500  E-O  003 
2.3015  6.56838806725i66E-obo3 
2.2989  5.999653198391796-0003 
2.2962  6i66533796303572E-0003 
2.2936  5.926929366963b6E-0003 
2. 2910  6f 1766bl6t699777E-0OO3 
2.2883  6.62026389001961E-0d03 
2.2857  6.217l6571178865E-0b03 
2.2831  6-.52258906595726E-0003 
2.2805  6.01287485275768E>0003 
2. 2779  6.26609328399117E-0003 
2.2753  6.32088590295865E-0003 
Z^nZl  6.66187687733222E-0003 
2.2701  6.617035203213571-0063 
2.2676  6^62952660170693E-0003 
2.2650  6. 29366116313568 E-0003 
2.2624  6.361929148852166-0003 
2.2599  6.165578119556872-0003 
2.2373  6.59406714178180E-0003 
2.2568  6. 33535296503851 E-0003 
2.2523  6.77716098312686E-0003 
2.2697  7. 0227091 1625273E-0003 
2. 2472  6.52612579373163E-0003 
2.2667  6.99302553012160E-0003 
2.2622  6.977696;32899997E-0003 
2.2396  6.95698055150523E-0003 
2. 2371  6.938038655597726-0003 
2.2366  6.52606198792311E-0003 
2.2321  6.81653OO99882816-O0o3 
2.2297  6. 79319926855015E-0003 
2.2272  7.04107091778639E-0003 
2.2267  6i66286765697067E-0b03 


2«2222 

2«2198 

ti2ni 

2.2148 
2c 2124 
2.2099 
2.2075 
2,Z05i 
2^2026 
2.2002 
2.1978 
2.1954 
2.1930 
2.1906 
2^1882 
2.1858 
2.1834 
2-1810 
2.1786 
2.1763 
2.1739 
2.1716 
2.1692 
2.1668 
2^1645 
2^1622 
2.1598 
2il575 
2il552 
2il529 
2*1505 
2.1482 
2.1459 
2.1436 
2.1413 
2.1390 
2.1368 
2.1345 
2.1322 
2.1299 
2.1277 
2.1254 
r.l23i 
2*1209 
2.1186 
2*1164 
2.1142 
2.1119 
2^1097 
2.1075 
2*1053 
2.1030 
2*1008 
2.0986 
2.0964 
2.0942 
2.0921 
1.0899 
2.0877’ 
2.0855 


6.666861627458336-0003 
6. 94 570 8 597 1 2219 E- ab  0 3 
7.0  023  8  8  549  0897  45-6003 

6.  9  999  2  0  76:9  2  24  8  4E-0  0  0  3 

7.  34801001561 181£-6o03 
7. 08  559  2103  0  599  2  E- 0  0  0  3 
7.  22377 127353 0l3E-b003 
6. 99229696634518E-6  003 
7.  71 43  0  032  03  5  8  8  0  E- 6  0  0  3 
T.  36  0 136  0 36 3 8358 E-0  0 0 3 
7.  27 3  5 0 1 3  0  9  7  8 97  6  E  - 6  0  0  3 
7.448335684514:97E-0  003 
7*  0418  5 1 63 3 8 43 5 2 E-O  0  0 3 
6. 84834662217781E-0003 
6.99779854329563E-0003 
7.  405468 177  6  232  2  E-O  0  0  3 
7.90655953923647E-0003 
7. 3  2746  5  916745275-0  0  0  3 
7.46673777203455E-0003 
7.86195253476762E-0003 
7^.  56  69  7  370  3  6  605  5  E-O  003 
7.54870525062046E-0003 
6.92096771538675E-0003 
7.  0935 19  0  31 8  6217  E-O  0  03 
6.6782006580297lE^b003 
6.83115982979388E-6o03 
7.064296982846945-0003 
7. 68  03759234991 3  E-OO  03 
7.75363184315836E-0003 
7.870ill4bl67238E^6003 
7.36036172898480E-0003 
7.60327272567451E-0003 
7.4629856996821BE^0003 
7. 0 9  5771746  8  7  5  8 1 E-O  0  0  3 
8.06644923874905E-0003 
7.56001017900587E-0003 
8.29679073810041E-0003 
7.63352866097478E-0003 
7.55352557573730Er0003 
8.23197375015638E-0003> 
8.15357882184742E-.0003 
fi  95987019891697E-0003> 
7.93043580688391E-0003 
7^669536289334648-000 3 
7.7925602'8284436E-:0003 
7.60021093051932E-0003> 
7.90565608048i55E*^0003 
8.05862470835450E-0003 
8.0694703839299lE-^0003 
8. 07873997074623 E-0003 
7.44014505334434Er0003 
7.845970422437i8E-0003 
7.820071829087285^0003 
8.32i60356522l38E-0003 
8 . 2 2793  6 10 1 2 63  8  2  E  ^0  0  0  3 
8.O561496B317059E-OOO3 
7.82771660489345E^0003 
8.48649372430637E^0003 
7.90810766304989E-0003 
8.09i86863e26028E^0003 


2.0833  8. 141^41 069823131-0003 

2.0812  7.9907  0711.0  332518-0003 

2.0790  7^.910573606025658-0003 

2.0768  8.282996279532998-0003 

2. 0747  8. 163671777978728-0003 

2.0725  !• 4639  20868 154628-0003 

2 . 07  0  4  7.5  04  5  82  0757  3246  8-00  03 

2.0683  .  7:. 853084^723066678-0 003 

2.0661  7.301750704108668-0003 

2 .0640  7.653813285  244  8 48- 000  3 

2 .0619  1^.63308481 0666768-0003 

2.0597  Z. 3  9 1 5  4  0  5 1 37370  38-0  003 
2.0576  71.295126193852758-0003 

2 .05  55  7.  3  3  0  8 1247  3 102198-00  0  3 

2.0534  6.90257638346026E-0603 

2.05 13  7 . 5  3 1 3  7 7  8  5  5 1 0  91 6  8-00  0 3 

2.0492  71.866635252895558-0003 

2 .0471  8 - 494 1749  5 1749768-00  0  3 

2.0450  8.637343545004238-0003 

2.042  9  8 . 468690 1 0 1935298-000  3 

2.0408  8.324239982499658-0603 

2.0387  8.303143431618378-0003 

2.0367  8.266423OO459923E-O0O3 

2.0346  8 . 8 169  877 1 5  39163  8-0  003 

2.0325  8.487066771337488-0603 

2.0305  9.208229234204178-0603 

2.0284  8^971635904373448-0003 

2.0263  8.679319777300528-0003 

2.0243  7i998504S0458450e-0003 

2.6222  8.160618163321968-0003 

2.0202  7.743142128482828-0003 

2.0182  6.861675514144858-0003 

2. 6161  6.685624206305368-0003 

2.0141  6.3 1084  5727114168-0  003 

2.0121  6.309051562006128-0603 

2.0101  6.889069909277628-0003 

2.0086  7.198391i9535298?8-6603 

2.0066  6.506601377624096-0603 

2.0040  5.733884252336458-0003 

2.0020  6.141754661278758-0003 

2-0000  6-278751463582928-0003 

l.i9980  7.952625420330iM-0003 

1-9960  8-603700755259076-0003 

1^9940  8^377988554713066-0003 

1.9920  8.394492058812386-0003 

1.9900  9.249123819088828-0003 

1.9881  8.752903749680746-0003 

1^9861  8.281826954451298-0603 

1.9841  8.409848598020638-0003 

1.9822  8. 855189988636886-0003 

1.9802  8.936925039634506-0003 

1^9782  8.475965058480038-0003 

1^9763  8.313105063308518-0003 

1^9741  8. 378974931048338-0003 

1.9724  8.474827493699888-6003 

149704  8.751266363177028-0003 

1.9685  8.746661093766768-0063 

149666  8. 1^2086753878178-6003 

149645  7.148438764836846-0003 

159627  7.418985445788458-6063 


1.9  6  08  8  •  6  0  221 6011 624 3 6 E*? 00  0 3 
l:.i9589  8>88t5032Tl08ll26^00(li 

1 .19563  8. 0825609i956i39E^0003 
t.^550  8.459696113803266-0  003 
1.9531  8\.  071515949548276^0003 
1.9  5 12  9. 18  5  3  0439  0  5  98  5  0  6^000  3 
1 .9493  9. 268161218429336-0003 

1.9474  .  8o937l514i629132E-0003 
1.94  55  9.2  189  17892392936-0003 

1.9436  8.777218333821066^0003 

1.94 17  9 . 61 099492  05 125  06-0003 

1.9399  9.908635475127886^0003 

1. 9380  9.857230734681376-0003 

1.9361  8. 978205296187986^0003 

1 . 9342  1.  0  5  88  8 19 1 7 1 33 196-0  002 
1.^9324  9. 865471327273676-0003 

1 .930  5  9.498551 1466063  6  6^0  0  03 
1^ 928  6  9. 48  73 14  9718  06  8  86-0  0  03 
1^9268  9.893640965726426^0003 


Appendix  8 
Figures  30-65 


93uaiaH!P  (uasjadi 


^inicrons 


100 


microns 


AN  Percent  Diffemce 


100- 


micrbhs 


8pu0ja|)!p  luasjsci 


microns 


microns 


microns 


MODTRAN  Percent  Difference 
Run  37 


r  <0 


89uajaii!p  tuaajad 


microns 


MODTRAN 


microns 


MODTRAN  Percent  Difference 
Run  06 


99U0J«)||p  |U99J9d 


microns 


MODTRAN  Percent  Difference 


r  (D 


I-  ■  I  I  • 


•ouaja^iip  luaaiad 


microns 


MOPTRAN  Percent  pifff»rence 
Run  14 


•ouaj9|j|!p  lusojad 


mlcronai 


MODTRAN  Percent  Difference 
Run  20 


ooudJ9)i|p  }uapjad: 


<ihicronsf 


^00  T  Ili  I  II  II mill  hi  jiit  MODTRAN  Percent  DIffe 


eouajaiMP  vudOJdd 


microns 


Figure  45 


microns 


M0btRAN  Percent  Dlfferefice 
Run  32 


f— 1— r— 1 — — I  r^r  — i--!  r  r'-i' 

8  g  S  §  g  °  ^  ^  8,  g  I 


03U9JOIIIP  ;}U89J8d 


ihiicronsj 


MODTRAN  Percent  Difference 
Run  36 


r  <D 


aouajaiiip  iuaojad 


sub'jaiui' 


100 


0OU8J8|){p  r)U80i8d 


90uajdii!p  luaojad 


microns 


•  <D 


•puaidiiip  )U89J8d 


4 


microns 


Figure  52 

LOWTRAN  Percent  Difference 


•pudJdiHP  Vuaojad 


microns 


microns 


tf> 


in 

« 

h 

9 

O' 

•H 

U4 


I — '  I  '•  r 


I  t  — r ”'“i — 

O  O  O  o  o 

^  8  «5*  T 


1  -r^  r-T-i  -I  f 


•  #  •  »  • 


/ 


eouajaiifp  luaajad 


microns 


83uajaii!P  )uaajad 


microns 


100  Tl! 


, microns 


iiriicrons 


opudjeilip  vUaojad 


Figure  59 


microns 


LQWTRAN  Percent  Difference 
Run  14 


microns 


Fi^re  61 

LQWTRAN  Percent  Difference 
Run  20 


83U8J9|}{p  luaojed 


microns 


o 


^  I  I — '  I  •  I 

^  ^  s  s 


8 


8 


T" 

8 


T-rr- 

8 


83U9jaji!P  luaojod 


microns 


LOWTRAN  Percent  Difference 
Run  32 


microns 


'Figure  65 

LOWTRAN  Percent  Piffefence 


microns 


